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Summary 

This manual describes the use of and relevant equations programmed in a computer code designed 
to carry out a comprehensive linear analysis of multilayered fiber composites. The analysis contains 
the essential features required to effectively design structural components made from fiber 
composites. The program is an outgrowth of two in-house computer codes, MFCA (Multilayered 
Filamentary Composite Analysis) and INHYD (Intraply Hybrid Composite Design). The inputs to 
the code are constituent material properties, factors reflecting the fabrication process, and composite 
geometry. The code performs micromechanics, macromechanics, and laminate analysis, including 
the hygrothermal response of fiber composites. The code outputs are the various ply and composite 
properties, composite structural response, and composite stress analysis results with details on 
failure. The code is in Fortran IV and can be used efficiently as a package in complex structural 
analysis programs. T he input-output format is described extensively through the use of a sample 
problem. The code manual consists of two parts. The mechanics for using the code are described in 
the first part, the pertinent equations programmed in the code are described in the second part. 


Introduction 

The importance of and need for a multilevel analysis used for designing structural components 
made of multilayered fiber composites are documented in reference 1. A multilevel analysis, which 
was efficient in predicting the structural response of multilayered fiber composites (with the 
constituent material properties, the fabrication process, and the composite geometry known), is also 
documented in reference 1. 

The multilayered analysis presented in reference 1 consists of (1) micromechanical theories for the 
thermoelastic properties and the stress-level limit of the single ply as functions of constituent material 
properties and the particular fabrication process, (2) the combined stress-strength criterion for the 
single ply, and (3) multilayered composite structural response and analysis (macromechanical or 
laminate analyses), where the interply layer effects are taken into account. A computer code designed 
to carry out this multilevel analysis, supplemented as noted by references 2 to 10, has been developed 
at the Lewis Research Center. This code is identified as MFCA for Multilayered Filamentary 
Composite Analysis (ref. 11). 

Intraply hybrid composites are a logical sequel to conventional and interply hybrid composites. 
Recently, theoretical and experimental investigations have been conducted on the mechanical 
behavior of intraply hybrids at the Lewis Research Center (refs. 12 to 14). The theoretical methods 
and equations described in these references, together with those for hygrothermal effects (ref. 15), 
have been integrated into a computer code for predicting hygral, thermal, and mechanical properties 
of intraply hybrid composites. This information can then be used in designing these composites. This 
code is identified as INHYD for Intraply Hybrid Composite Design (refs. 16 and 17). 

The present computer code is a synergistic combination of the aforementioned computer programs 
MFCA and INHYD together with several significant enhancements. The code is referred to as ICAN 
for Integrated Composite Analyzer. It utilizes the micromechanical design of INHYD and the 
laminate analysis of MFCA to build a comprehensive analysis capability for structural composites. 
Additional features unique to ICAN are the following: 



(1) Ply stress-strain influence coefficients 

(2) Microstresses and microstress influence coefficients 

(3) Stress concentration factors around a circular hole 

(4) Calculation of probable delamination locations around a circular hole 

(5) Poisson’s ratio mismatch details near a straight free edge 

(6) Free-edge stresses 

(7) Material cards for finite-element analysis using NASTRAN or MARC 

(8) Failure loads, summary based on the maximum stress criterion and laminate failure stresses, 
and summary based on first-ply failure and fiber breakage criteria 

(9) Transverse shear stresses and normal stresses 

In addition to the above, ICAN has its own data base of material properties for commonly used 
fibers and matrices. The user needs to specify only the coded names for the constituents. The 
program searches and selects the appropriate properties from its library. Furthermore, input data 
preparation has been simplified substantially by the introduction of a partial free-field format. The 
output formats have also been improved significantly to ease user interpretation of the results. These 
enhancements make ICAN significantly more user friendly than its predecessors. The computer code 
has been programmed in Fortran IV and has been tested in UNIVAC 1108, IBM 370, and CRAY 1 
computers. 

Since this report is to serve as a users manual, the code is divided into two parts, the users manual 
and the programmers manual. The Users Manual describes the mechanics of using the code with 
respect to program format, input and output, and sample input data sets. The descriptions are 
extensive enough sc that even designers and analysts with little or no programming experience can 
easily use the code. 

The programmers manual gives the various subroutine descriptions and the equations programmed 
therein, with details on the input and output and the global storage locations. This, along with the 
listing of the source program, allows the user to make his own modifications to the code as they 
become appropriate for further enhancements. 

The Fortran variables are defined in appendix A. Included is information such as which part of the 
program of each global variable is generated. Table I provides a summary of details for preparing 
data cards, and the input data given in table II provide for immediate testing of the code. Properties 
for a few commonly used fibers and matrix materials are listed in appendix B. Appendix C shows 
sample input and output data for a specific case. 


Symbols 


A cx 

composite axial stiffness 

aR 

74 or 

reduced axial stiffness 

BIDE 

boolean, true if interply effects are included 

Ccx 

composite coupling stiffness 

C e i 

string with force variables 

C e 2 

string with displacement variables 

COMSAT 

boolean, true if laminate analysis is wanted 

CSANB 

boolean, true if membrane and axial symmetry exists 

D C ,D( 

moisture diffusivity 

Dcx 

composite flexural rigidity 

dr 

‘-'cx 

reduced bending rigidity 


displacement vector 

d f 

filament equivalent diameter 

Ef , E c f 

filament elastic constants 

Efu, etc. 

fiber normal modulus 

E ( ,E C( 

ply elastic constants 


E n{ , etc. 
Em* E C m 
Efn \ 1 » 
&m,etc 

F 

Gyi2, etc. 
Gm» etc - 

G m\2 > etc * 

H, 


H 


kc 


u 


K, 


cl 1 ,c22,c33 


K t 


cxx, cyy, cxy 


A /ll 

Km 11 

k v 

k f 

kf ( 


k ml 

k v( 

k sc 

M ( 
\f rY 


M, 


cMsX 


M. 


cT f x 

m 

Hex 

^ cMfX 
^cTfX 

N f 

N ( 

Nk 

NONUDF 


N, 


pc 


N, 


P? 


1 cp 

Pi 

Pip 

Qf,i,p t r,s 

R 

RINDV 


ply normal modulus 

matrix elastic constants 

matrix normal modulus 

matrix failure strain allowables 

combined stress-failure function 

fiber shear modulus 

ply shear modulus 

matrix shear modulus 

interply distortion energy coefficient 

array of constituent heat conductivities 

composite heat capacity 

index, generally ply or interply 

composite three-dimensional heat conductivities 

composite two-dimensional heat conductivities 

fiber heat conductivity 

ply heat conductivity 

matrix heat conductivity 

apparent void volume ratio 

actual fiber volume ratio 

ply apparent fiber volume ratio 

actual matrix volume ratio 

ply apparent matrix volume ratio 

ply apparent void volume ratio 

array of limiting conditions 

ply moisture 

applied moment 

hygral moment 

thermal moment 

load condition index 

applied membrane loads 

hygral force 

thermal force 

number of fibers per end 

number of plies 

number of load conditions 

boolean; true if detailed Poisson’s ratio differences chart is to be suppressed 

string PROPC length 

string PROP length 

composite properties array 

string PROPC 

ply properties array 

string PROP main program 

indices to print out string PROP 

transformation matrix 

boolean; true if displacements are known 


s c 

S n i7-, etc. 

h 

ft 

w cb 

x,y,z 

<*c 

«/ 

a ( 

a m 

Pc 

P e >Pc 

Ph 

Pt>Pf>@m 

Ps 

Pv 

5 / 

6 f 

5, 

e cs 

6 CSX 

<T 

Qti’Qfc 
vf\i, etc. 
v n2 , etc. 

v m\2' etc. 

Pf.m 

Pmw 

1,2,3 


composite failure stress 
ply limit failure stresses 
ply temperature 
ply thickness 

composite local curvature changes 

structural reference axes 

composite coefficient of thermal expansion 

fiber thermal coefficient of expansion 

ply thermal coefficient of expansion 

matrix thermal coefficient of expansion 

moisture expansion coefficients 

correlation factors for ply thermoelastic properties 

correlation factor for ply heat conductivity 

moisture expansion coefficients for ply, fiber, and matrix 

correlation factor for ply strength 

matrix strain magnification due to ply strain in the presence of voids 
interfiber spacing 
interply layer thickness 
interfiber spacing 

angle between composite material and structural axes 
reference plane membrane strain 
ply strain 

angle between ply material and composite axes 

fiber Poisson’s ratio 

ply Poisson’s ratio 

matrix Poisson’s ratio 

fiber and matrix weight density 

density of matrix with moisture 

ply stresses, fiber stresses, and matrix stresses 

material reference axes 


Users Manual 

The mechanics required to use this code for the analysis of multilayered fiber composites are 
described in this part of the report. The theory on which the code is based is described in the second 
part of the report (Programmers Manual). 

The physical representations of the constituents used in the code are illustrated in figure 1 . This 
figure shows a complete integration schematic starting with the constituent materials, fiber and 
matrix. The required input properties and computed properties at various levels are summarized in 
symbolic form as follows: 

( 1 ) Properties required by code as input for a fiber: Ef\\ f 22 33 » ^f\2 23 13 > ^-*f\2 22 13* ®-f\ \ 2233 * 

^/*l 1,22,33; Hcf* Pf ; Nf \ and S/r 

(2) Properties required by code as input for a matrix: £^ 11 ,22,33* ^ml2, 23,13* ^/12,23,13 a /wl 1,22,33* 

11,22,33* Pnv ® mpt * &mpc* ^ mps* an d £>mpTOR' 

(3) Properties required by code as input for a single ply: fiber and matrix properties and ply 
characteristics f3 m fly, and 7). 




Figure 1. — Schematic of typical multilayered fiber composite and some basic components. 


(4) Properties computed by code for single ply: £,’11,22,33; ^12,23,135 ^12,23,13; 1,22,335 22, 335 

H c f Pri 6?; 11c, 227',22C,i25, 2355 K(\2> anc * stress analysis factors 6^1,22,12; a n 1,22,12; an< ^ 

1.0 -F(<j,S,K m ). 

(5) Properties required by code as input for a composite: ply properties and composite 
characteristics Hj, K? [2a p, N cv , M cx or U cx , and W cx . 

(6) Output computed by code for a composite: [e cx \ = [£* c ]{a c J + 7Via c *j; l^cl / Kcxx,vy,xy> H c \ 


Hex 


A CX^CX 


Ucx 

! 

N cxT e 





1 

+ 

1 


M cx 


C D 

y “'CX L 'CX 


W 1 

fY cx t 

1 1 

M C xT, 


and the inverse A<pd e ij* 

Figure 2(a) shows the subroutines and sequence in the code. The subroutines between the Main 
program and the input data may be arranged in any desired order. The user should refer to figure 2(b) 
for the logic flow of the analysis. 

The following four steps are required to use the code in the user’s computer facility: 

(1) Obtain the code 

(2) Make it operational in the user’s computer facility 

(3) Supply the input data 

(4) Interpret the code output results 

Obtain the Code 

The code may be obtained in cards. If this is not convenient or possible, the cards can be punched 
from the compiled listing (contact COSMIC, The University of Georgia, Athens, GA 30602, 
concerning the availability of this program). 
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JCL cards 



Main program, ICAN 

2 



Micromechanics, INHYD 



Composite properties, FIBMT 

4 



Hygrothermomechanical 
response, HTM 

|5 


Hybrid composite properties, 
COMPP 



Maximum of an array, AMAXF 


Three-dimensional composite 

7 

Failure mode and failure load. 

20 

properties, GACD3 


FLBLD 




8 


21 

Two-dimensional composite 

Microstress influence 

properties, GPCFD2 


coefficients, MINCOF 



Flexural strengths, FLEXX 

9 

Name strings, BLOCK DATA 

22 

Composite stress analysis, 
COMSA 

10 

Print the emblem ICAN, LOGO 

" 23 



11 


24 

Read constituent properties, 


User input, Input data 


BANKRD 


unit 5 



Calculate free edge stresses, 

12 

Data base, FBMTDATA. BANK 

EDGSTR 


unit 8 


Property generator for 
INHYD, IDGER 


(a) Schematic showing subroutines and sequence of ICAN code. 
Figure 2. — Subroutines, sequence, and logic flaw of ICAN code. 


Make it Operational 

A prerequisite to the program is the availability of a Fortran compiler in the user's computer 
facility. To run the program, certain computer-system-dependent control cards (Job Control 
Language (JCL) cards) may also be necessary. The computer system personnel should be consulted 
about these items. 

Once the deck of cards has been assembled (except input data) with the proper control cards as 
shown in figure 2, the user is ready to compile the code in his facility. The compilation will indicate 
whether any additional modifications are needed. Most modifications will be minor and will usually 
deal with certain Fortran statements peculiar to each compiler. 

Supply the Input Data 

The physical arrangement of the input data cards is illustrated in figure 3. Details for preparing the 
input data cards are summarized in table I. A detailed description of these cards is given 
subsequently. A sample for preparing input data for a four-ply symmetric laminate is presented in 
table II. This laminate has two different material systems. The 0° plies are of AS graphite 
fiber/intermediate-modulus, low-strength epoxy matrix composite. The 90° plies are made of a 
hybrid composite. The primary composite is S glass/high-modulus, high-strength epoxy. The 
secondary composite is AS graphite/intermediate-modulus, high-strength epoxy. 






(b) Schematic showing logic flow of ICAN code. 
Figure 2. — Concluded. 



Figure 3.— Physical arrangement of input data cards. 
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TABLE I. -SUMMARY OF DETAILS FOR PREPARING INPUT DATA CARDS 


Card 

group 

Identification 

Code symbol 

Number of 
entries 

List of 
entries, 
sequential 
order 

Card 

field 

columns 

Format 

Comments and 
engineering units 

1 

Title card 

TITLE 

80 

Alphabetic 

characters 

1 to 80 

(10a8) 


2 

STDATA 

NL.NLC.NMS 

3 

N t .N t ,,N, m 

9 to 36 

(a8,3I8) 

Composite geometries 

3 

Boolean for input 
displacement 
Boolean for interply 
layer energy 
contribution 
Boolean for mem- 
brane and bending 
Symmetry- 

Boolean for laminate 
analysis 

Boolean for Poisson 
ratio chart 

RINDV 

BIDE 

CSANB 

COMSAT 

NONUDF 

1 

1 

I 

I 

1 


1 to 6 
1 to 6 

1 to 6 

1 to 6 
1 to 6 


T (true) if displays are inputs; 
otherwise F (false) 

T (true) if contributions are 
desired; otherwise F (false) 

T (true) if symmetry exists; 
otherwise F (false) 

T (true) if laminate analysis is 
desired; otherwise F (false) 

T (true) if Poisson ratio differ- 
ences chart is not desired; 
otherwise F (false) 




(L6) 






4 

PLY (ply desired) 

INPl.IPl.TU.TCU 

DELM.TETA.THCKNS 

7 

U'T*, 

T CU ,AM, 

0 ( ,tf 

1 to 64 

(a8.2I8,5F8.3) 

Ply layup and temperature 
and moisture conditions 

5 

MATCRD (material 
system details) 

CODES(LLJ), 

CODES(l,2,J), 

VFP.VVP, 

CODES(2,LJ). 

CODES(2,2,J), 

VFS.VVS.VSC 

9 

Primary 

composite 

code 

names for 
fiber and 
matrix, kf, 

k v ; 

Secondary 
composite 
code 
names 
for fiber 
and 

matrix, 

k.sc> kft k v 

I to 64 

(a8,2a4,2E8.2) 

Description of material 
systems to be used 

6 

PLOAD (loading 
details) 

NX,NY,NXY,THCS 

MX.MY.MXY 

i 

DMX.DMY 

PRSSU.PRSSL 

4 

3 

4 

N X ,N,,N XV , 

0 C sM x .M„ 

\M„ 
dM x /</,. 
dM y ld y . 

p«A 

1 to 32 
1 to 32 

1 to 40 

(a8,7E8.4) 

(a8,7E8.4) 

(a8,7E8.4) 

Loading conditions (inplane) 
Angle of inclination to the 
structural jt-axis 
Loading conditions (bending) 
Loading conditions 
(tranverse) 


Input data for additional composite systems may be easily prepared. This is done by selecting the 
desired fiber and matrix from the available materials listed in appendix C using the variable 
FBMTDATA.BANK and modifying the appropriate entries in the input data sample illustration. 

After the input data have been properly assembled (as shown in fig. 3), they are placed in their 
physical position (fig. 2) and the code is ready to be run. 

Detailed Description of Input Data 

The card group numbers referred to here are given in figure 3 and table I. The sequential order of 
the entries in each card group is given in table I. Note that most data cards are identified by a 
mnemonic to indicate the card group in which it belongs in the input data deck. Also, most datacards 
are divided into fields of eight, with one entry per field being allowed. The mnemonic is entered in 
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TABLE II.-1CAN SAMPLE INPUT DATASET 


FOUR PLY SYMMETRIC LAMINATE 


ICAN SAMPLE INPUT 

DATA. 



STDATA 

4 

1 


2 




T 





COMSAT 



F 





CSANB 



F 





BIDE 



F 





RINDV 



T 





NONUDF 



PLY 

1 

1 

70.00 

70.0 

0.0 

0.0 

0.010 

PLY 

2 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

3 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

4 

1 

70.00 

70.0 

.0 

.0 

.010 

MATCRDAS- - IMLS 

0.55 


0.02 

AS-- IMLS 

0.0 

0.57 

0.03 

MATCRDSGLAHMHS 

.55 


.01 

AS--IMHS 

.4 

.57 

.01 

PLOAD 1000. 

0.0 


0.0 

0.0 


NX , NY , NXY , THCS 

PL0AD 0.0 

.0 


.0 



MX, MY 

MXY 

PLOAD .0 

.0 





DMX/QX , DMY/QY , PRSS 


format A8, and the integers are entered in format 18. The real numbers may be entered anywhere in 
the appropriate field. The following is a brief description of each card group together with examples 
taken from table II: 

Title card . 


Any title of length up to 80 characters 


FOUR PLY SYMMETRIC LAMINATE. ICAN SAMPLE INPUT DATA. 


As shown, any title of length up to 80 characters including blanks may be supplied on this card. 
Starting data card . 


I 8.9 16,17 24,25 a 32 


Mnemonic 

N t 

A'V 

N fm 

STDATA 

4 

] 

7 


This card has a mnemonic STDATA. It contains the overall laminate and loading details. These 
details are the number of plies N$, the number of loading conditions and the number ot different 
material systems N ms . 

Boo leans. 


I 6.7 


Boolean 
T or F 

This space may be used for comments 

T 

COMSAT 

F 

RINDV 

F 

BIDE 

F 

CSANB 

T 

NONUDF 


A set of booleans, COMSAT, RINDV, BIDE, CSANB, and NONUDF is defined through these 
cards. These are five cards, one per each logical variable. The format is L6. The variables have the 
following functions: 

(a) COMSAT.— The letter T in the card will direct the program to perform a complete laminate 
analysis. A letter F would terminate the program at this stage. 







(b) RINDV.— The letter T is entered in the card if the displacements are inputs; otherwise, the 
letter F is entered. 

(c) BIDE. — The letter T is entered in the card if the interply layer contributions on the composite 
are desired; otherwise, the letter F is entered. 

(d) CSANB. — The letter T is entered in the card if the composite has both membrane and bending 
symmetry; otherwise, the letter F is entered. 

(e) NONUDF. — The letter T is entered if the detailed Poisson’s ratio difference chart is to be 
suppressed; otherwise, the letter F is entered. 

Ply details card group. 


1 8,9 16,17 24.25 32.33 40,41 48,49 56,57 64 


Mnemonic 

Ply 

Material 

MID 

T it 

T 

1 cu 

M 

e, 

U 

PLY 

1 

1 

70.00 

70.0 

.0 

0.0 

.015 

PLY 

-> 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

3 

2 

70.00 

70.0 

.0 

90.0 

.005 

PLY 

4 

1 

70.00 

70.0 

.0 

0.0 

.010 


All the cards in this group have the mnemonic PLY. The number of cards is with eight entries on 
each card. The first entry is PLY. The second and third entries are identification numbers for the ply 
and the material system, respectively. The fourth and fifth entries are the use temperature T u and the 
cure temperature T cu , respectively. The sixth entry is the percentage of moisture M. The seventh and 
the eighth entries are the orientation angle 6 of the ply and the thickness of the ply, respectively. A 
default value of 0.005 is taken for the thickness if this entry is missing. The material system 
identification number should be different not only for different composite systems but also for 
varying use temperature or moisture content from ply to ply. 

Material system details. 


1 8,9 16,17 24,25 32,33 40,41 48,49 56,57 64 


Mnemonic 

Fiber, 

matrix 

k f 

k v 

Fiber, 

matrix 

n<- 

k r 


MATCRD 

AS— IMLS 

.55 | 

.02 

AS— IMLS 

0.0 

.57 

.03 

MATCRD 

SGLAHMHS 

.55 

.01 

AS— IMHS 

0.4 

.57 

.01 


All the cards in this group have the mnemonic MATCRD. The number of cards is N ms with 10 entries 
in each card. The first entry is MATCRD. The second and the third entries are coded words for fiber 
and matrix material of the primary composite. The code words are entered in format 2A4. For 
example, the code for AS-type fiber is AS— and epoxy matrix is EPOX. A dictionary of codes for 
several fibers and matrices is provided in appendix C. The user may choose any combination of fiber 
and matrix for a composite system. The fourth and the fifth entries pertain to the details of the 
primary composite system. They are the primary fiber volume ratio and the primary void volume 
ratio, respectively. The next two entries refer to the secondary composite system which is applicable 
for the case of the hybrid composite ply. They should be the same as the second and third entries for 
standard composite systems. The next entry is the secondary composite system volume ratio. This is 
zero for the standard composite systems. The last two entries are the fiber volume ratio and the void 
volume ratio for the secondary composite system. These values are entered when applicable. 

Load cards . 


1 8,9 16,17 24,25 32,33 40 


Mnemonic 

N* 



T hts 

; PLOAD 

1000. 

0.0 

0.0 

0.0 

PLOAD | 

0.0 

0.0 

0.0 


PLOAD 

0.0 

0.0 

0.0 

0.0 





All the cards in this group start with the mnemonic PLOAD. There are three cards for each loading 
condition. Thus, the total number of cards is 3/Vfc.. The first card under each loading condition 
contains entries N x , 7V v , and N xy for the membrane loads and Thcs for the orientation of the loads 
with respect to the structural axes. Similarly the second card contains the bending resultants M x , A/ v , 
and M xv . The last card contains the transverse shear resultants DM X and DM y and the transverse 
pressures P u and P ( . 

The user input data are read from I/O unit 5. Apart from this, ICAN uses two more units, 7 and 8, 
for its I/O operations. Unit 8 is used to store the material properties data base. Unit 7 is used as a 
scratch file by ICAN. These I/O units must be appropriately defined by using the operating system 
JCL. 

Output 

The following items are printed out by the program: 

(1) ICAN logo 

(2) ICAN coordinate systems 

(3) ICAN input data echo 

(4) Input data summary 

(5) Fiber, matrix, and ply level properties of primary and secondary composites 

(6) Composite three-dimensional strain-stress and stress-strain relations about the structural axes; 
MAT9 card for MSC/NASTRAN solid elements 

(7) Composite properties generated in array PC 

(8) Composite constitutive equations about the structural axes 

(9) Reduced bending and axial stiffnesses 

(10) Data for finite-element analysis 

(11) Displacement-force relations for the current load condition 

(12) Ply hygrothermomechanical properties/response 

(13) Details of Poisson’s ratio mismatch among the plies 

(14) Free edge stresses 

(15) Microstresses and microstress influence coefficients for each different composite material 
system 

(16) Stress concentration factors around a circular hole 

(17) Locations of probable delamination around circular holes 

(18) Ply stress and strain influence coefficients 

(19) Laminate failure load analysis based on the first-ply failure/maximum stress criteria 

(20) Summary of the laminate failure stress analysis based on two alternatives, first-ply failure and 
fiber breakage 

The printout of the input data summary (app. B item 4) shows details regarding composite 
geometry, constituent specifications, temperature and moisture conditions, and the loading 
conditions. 

The next few pages of the output are generated by the INHYD program package. They show the 
fiber-matrix properties for the different composite systems and the ply level properties of the 
composites (app. B, item 5). 

The output of the composite three-dimensional strain-stress temperature and moisture relations 
and composite stress-strain relations about the structural axes are printed under the following 
headings: 

(a) 3-D COMPOSITE STRAIN STRESS TEMPERATURE MOISTURE RELATIONS- 
STRUCTURAL AXES 

The matrices [E c \g and [/3 c j 5 in the equation 

l 6 cl5 = [E c \s — ~~Mf[f} c ] s 

where A T ( = T ( - T cu 

are printed by the subroutine GACD3. 



(b) 3-D COMPOSITE STRESS STRAIN RELATIONS-STRUCTURAL AXES 
The matrix in the equation 

!°cls — l-^cJ.rl^cSj 

is printed out by the subroutine GACD3. 

The subscripts in the preceding equations indicate that the relations are written about the structural 
axes. It is noted that these properties are only local to subroutine GACD3. They can be made global 
if needed. The properties needed to prepare the MAT9 card of MSC/NASTRAN are printed out next 
under the heading MAT9 CARD FOR MSC/NASTRAN SOLID ELEMENTS (app. B, item 6). 

The output of the composite properties, generated in array PC, are printed under the following 
heading (app. B, item 7): 

COMPOSITE PROPERTIES— VALID ONLY FOR CONSTANT TEMPERATURE THROUGH 
THICKNESS LINES 1 to 31: 3-D COMPOSITE PROPERTIES ABOUT MATERIAL AXES 
LINES 33 to 62: 2-D COMPOSITE PROPERTIES ABOUT STRUCTURAL AXES 
Sixty-two entries are printed under this heading as shown in the following list: 


Code name 

Notation 

Explanation 

PC(1) 

Pc 

weight density 

PC(2) 

‘c 

thickness 

PC(3) to PC(ll) 

\£ c \ 

three-dimensional stress-strain relations about material axes 

PC(12) to PC(14) 


three-dimensional coefficients of expansion about material 



axes 

PC(15) to PC(18) 

\K C \,H C 

three-dimensional heat conductivity and heat capacity along 



material axes 

PC(19) to PC(30) 

Ec\bGc\2’ v i\2 

three-dimensional constants about material axes 

PC(31) 

Z c 

distance to reference plane from bottom of composite 

PC(32) 


blank 

PC(33) to PC(38) 

\E C \ ~ 1 

two-dimensional stress-strain relations about structural axes 

PC(39) to PC(47) 

Ec\\’ G c\2’ v c\2 

two-dimensional elastic constants along structural axes 

PC(48) to PC(54) 

a c , K c , H c 

two-dimensional coefficients of thermal expansion, heat 



conductivity, and heat capacity along structural axes 

PC(55) to PC(58) 

D c 

moisture dif fusi vities 

PC(59) to PC(62) 

0c 

moisture expansion coefficients 


Array PC and its corresponding string and headings are controlled by the formats in subroutine 
GPCFD2. For nonuniform temperature/moisture, the bending equivalent and the membrane 
equivalent elastic constants may be obtained by utilizing the reduced bending rigidity matrix and the 
reduced stiffness matrix which are also regular output of ICAN. 

The output for the composite constitutive equations are printed in the following manner (app. B, 
item 8): 


FORCES FORCE DISPLACEMENT DISPL T-FORCES H-FORCES 

RELATIONS 


i*«j 


[A cx \\C cx \ 



( \ N CT r x\ , 

i ; 

\ n cm ( x\ 

\M CX \ ' 

1 

\C cx \\D cx \ 


\w cb \ 

| ( \ m ct,x\ 

r ' 

1 1 

WcM t x\ 


The elements of matrices A cx , C ex , D cx , N cT N cM M cT and M cMx are printed out by the 
subroutine GPCFD2. 

The output for the reduced bending rigidities is printed under the heading (app. B, item 9): 
REDUCED BENDING RIGIDITIES. The elements of [£&] are printed out as a matrix. 

Similarly, the output for the reduced axial stiffness [A^ x ] is printed out under the heading 
REDUCED STIFFNESS MATRIX. The corresponding formats for the above two outputs are in 
subroutine GPCFD2 (app. B, item 10). 


The next printout comes from the main program under the heading: SOME USEFUL DATA FOR 
F.E. ANALYSIS. This information is useful for preparing material data cards for finite element 
codes NASTRAN and MARC. 

The inverse of the constitutive equations is printed out in the following manner (app. B, item 11): 


DISP DISPLACEMENT FORCE FORCES 

RELATIONS 


i^caxl 


1 — 

ft 

i 

1 

, l*«l + I*c7}*l 


K*l ] 


I Qvl \D C x\ 

i 

\{M CX \ \M cT(X \ 

\ M cM ( xl 


The elements of this inverse are printed out in the subroutine COMSA. 

The current values for the loads and corresponding set of ply properties generated in array PL are 
printed out next (app. B, item 12). The explanations of the 75 entries in the PL property array are 
given in the following list: 


Code 

name 

Notation 

Explanation 

PL(1,I) 


k v 

ply void volume ratio 

PL(2,I) 


kf( 

ply apparent fiber volume ratio 

PL(3,I) 


k f 

ply actual fiber volume ratio 

PL(4,I) 


k mS 

ply apparent matrix volume ratio 

PL(5,I) 


k m 

ply actual matrix volume ratio 

PL(6,I) 


Pt 

ply weight density 

PL(7,I) 



ply layer thickness 

PL(8,I) 



ply and interply layer thickness 

PL(9,I) 


Hj 

interply layer distortion energy coefficient 

PL( 10,1) 


Z( 

distance from bottom of composite to ply centroid 

PL(1 1,1) 


z cg 

distance from reference plane to ply centroid 

PL(12,I) 


Ocs 

angle from structural axes to composite material axes 
(same for all plies) (fig. 2) 

PL(13,I) 


0( 

angle from ply material axes to composite material axes 
(fig. 2) 

PL(14,I) 



angle from ply material axes to composite structural axes 
(fig. 2) 

PL( 15,1) to 

PL(23,I) 

[£(■] - 1 

ply stress-strain relations 

PL(24,I) to 

PL(26,I) 

l«fj 

ply thermal coefficients of expansion 

r 

73 

o 

PL(29,I) 

[X(\ 

ply heat conductivities 

PL(30,I) 


Hot 

ply heat capacity 

PL(34,I) to 

PL(42,I) 


ply elastic constants 

*0 

r 

OJ 

O 

PL(48,I) 

D ( and /3 ( > 

moisture diffusivities and expansion coefficients 

PL .(49, 1) 



interply delamination factor 

PL(50,I) 


T, 

ply temperature 

PL(51 ,1) to 

PL(60,I) 

^(’1 1 T' etc ' 

ply limiting stresses 

PL(61,1) 


K n2a& 

coefficient in combined stress-strength criterion 

PL(62,I) 


combined stress-strength criterion 

PL(63,I) 



interply delamination criterion 

PL(64,I) to 

PL(69,I) 

M.M 

ply applied strains and stresses 

PL(70,I) 


Apj 

adjacent ply relative rotation 

PL(71 ,1) 


Hoffman’s failure criterion 

PL(72,I) 



ply moisture 

PL(73,I) 


a m 

transverse shear stress 

PL(74,I) 


0(23 

transverse shear stress 

PL(75,I) 


a (33 

thickness stretch stress 



The next printout shows Poisson’s ratio differences between the plies and the composite (app. B, 
item 13). They are printed out by the subroutine FESTRE under the heading DETAILS OF 
POISSON’S RATIO MISMATCH. 

The stress peaks near the free edge region are printed out next by the subroutine EDGSTR under 
the heading (app. C, item 14) FREE EDGE STRESSES. 

Item 14 of appendix B shows ply stresses in the structural coordinate system and the through-the- 
thickness stresses a zz , a xz , and o yz . The boundary layer decay length is also shown in the table under 
the heading YDCAY LENGTH. Care must be exercised in interpreting the results. They are based on 
approximate engineering theories and give good qualitative information regarding the relative 
magnitudes of the peaks in the individual plies. This printout is suppressed in the case of combined 
loading. 

The microstresses in each ply are printed out next by the subroutine MCRSTR (app. B, item 15(a)). 
Two regions of interest are considered for the computations, the region between the fibers composed 
entirely of matrix (A) and the region consisting of fibers as well as matrix (B). The stresses are given a 
descriptive notation. Thus, SM2AL means stress in matrix along the transverse (2) direction in region 
A due to a ply stress along the longitudinal direction of the material. Figure 4 shows the definitions 
for regions A and B. The printout also shows microstresses resulting from moisture and temperature 
differences if nontrivial Afyand 7} are present. 

The microstress influence coefficients, stresses due to unit applied stresses in direction 11, 22, 12, 
13, and 23 (app. B, item 15(b)); unit temperature difference Tf y and unit moisture content A/^are 
output from the subroutine MINCOF. These variables are printed out following the microstresses. 

Under the heading STRESS CONCENTRATION FACTORS (app. B, item 16) are printed out the 
factors K\ yy , and xy which are due to inplane loading around a circular hole at 5° intervals by 
the subroutine STRCNF. Cumulative stress concentration due to combined loading may be estimated 
by simple addition of the respective stress concentration factors. 

The next output (app. B, item 17) is under the heading POISSON RATIO DIFFERENCES and 
results from the subroutine NUDIFS. For each ply, the Poisson’s ratio differences, (j^-pj. -1 ) and 
anc * the products K^ xx (v\— v c ) y K\ yy (tf- v c ), and K\ xy (*4 _ *v) are printed out at 0 intervals 
of 5° around a circular hole. This is suppressed if the boolean NONUDF is set to TRUE. This item 
shows the locations of probable delamination for each ply. These are the locations where products 
such as K\ xx 04“ y/)> for example, are maximum. 

r Pty(£) 



Figure 4. — Definitions of regions for ply microstress calculations. 



The next item in appendix B shows the ply stress and strain influence coefficient arrays and ply 
stress influence coefficient arrays (app. B, item 18). These are computed in the subroutine COMSA 
and are printed by the main program ICAN. The first table gives the influence coefficients based on 
unit loads or moments/inch. The second table gives the influence coefficients in terms of unit applied 
stresses. Explanations of usage of these tables are provided at the end of each table. 

The output from the subroutine MSCBFL is printed out next under the heading LAMINATE 
FAILURE STRESS ANALYSIS (app. B, item 19). The analysis is based on first-ply failure criteria. 
Results are printed in a tabular form for each ply, and a summary of the analysis is shown in the end 
(app. B, item 20). The summary shows the critical ply, the failure mode, and the load for each of the 
applied load types, o cxxT , o cxxC , o cyyT , o cyyCi and o cxy s „ respectively. The first table shows results 
based on first-ply failure, and the second table shows results based on fiber failure by breakage. 


A Typical IBM Terminal Session 

To run ICAN, the user must first install and compile the program on his/her computer according 
to the system to be used. The procedure used on the Lewis Research Center IBM 370 is described in 
detail here starting from log on. The computer prompt signals are identified with uppercase letters. 
User entries are in lowercase letters. The following are prerequisites for the user to be able to run 
ICAN: 

(1) A knowledge of how to compile and store the object module in the public storage space. 

(2) A knowledge of redit or tedit processors so as to be able to create vs datasets of the input data 
deck. The details of the input data format have already been described in earlier paragraphs. 

(3) A knowledge of commands like rmds, mds, ddef, libdef, print, and erase. These are a few of the 
commands commonly used in running a program on the IBM 370. 

The user is advised to migrate the object deck, the input dataset, and the material property data 
base so as to conserve his/her permanent storage. The object deck, which is a binary version of the 
compiled source program, is referred to here as OBJ. ICAN. The data base of material properties is 
referred to as FBMTDATA.BANK. 

The session is started by logging on at the terminal. This is achieved by typing logon, userid, and 
password. The system replies 

TSS/370 RELEASE 3.0 PRPZ3 FTF18 

SOME MESSAGE TASKID = OBD7 POOLID = LRCFM - 

LOGON AT 11:30 ON 01/15/84 

The user is now ready for the session. The first phase of the session consists of restoring the 
necessary data sets to temporary storage. This is achieved by the following commands: 

rmds obj. ican, aaa 

SUCCESSFUL (TEMP) RESTORE OBJ. ICAN AS (AAA) 
rmds fbmtdata.bank, ccc 

SUCCESSFUL (TEMP) RESTORE FBMTDATA.BANK AS (CCC) 
rmds ican. sample. input, bbb 

SUCCESSFUL (TEMP) RESTORE ICAN. SAMPLE. INPUT AS (BBB) 

At this point, the user has all the necessary data sets to run ICAN in his/her temporary storage. 

The input/output fortran units that are utilized by ICAN for its various input/output operations 
need to be defined next. This forms the second phase of the session and is achieved by 

ddef ft05f001, vs, bbb 
ddef ft06f001, vs, icanout.bbb, ret = t 
ddef ft08f001, vs, ccc, ret = t 
ddef ft07f001, vs, T7, ret = t 

During these operations, the system usually responds by the minimum prompt, the underscore (_). 

The third phase consists of loading and executing the object deck and printing out the results. This 
is done by the following commands: 



libdef Ids, aaa 
load gpcom$$$ 
ican 

TERMINATED: STOP 
print icanout,bbb, prtsp=edit 
PRINT BSN = 8835, 1200 LINES 

The last phase involves cleaning up the user’s storage place and is achieved by issuing the 
commands 

release Ids 
release ft 
erase aaa 
erase bbb 
erase ccc 
erase t7 

The user may either logoff or proceed to execute another run for a different set of input data after the 
preceding set of commands. 


Programmers Manual 

A brief description of the main program (or control program) and theoretical equations 
programmed in the code are presented in this portion of the report. The subroutine descriptions 
follow the order of execution as shown in the flowchart (fig. 2(b)) rather than the physical sequential 
order (fig. 2(a)). It is assumed in the following discussion that the user has a working knowledge of 
computer programming and that he/she is familiar with the terminology appropriate to multilayered 
composite mechanics. 

The assumptions and details leading to the derivation of the equations programmed in the code are 
not included here. However, they are described in the references cited. It is suggested that the 
interested user have these references available to him/her. 

The information provided in this portion of the code together with the source program listing 
enables the user to modify, implement, and extend the code according to need. 


Main Program 

The main program contains the global variables, the various subroutines, the input data and 
format, the various program control statements, and the output. These are discussed subsequently. 
The flowchart of the program is shown in figure 5. 

The global variables are given in the following list: 

boolean CSANB, BIDE, RINDV, COMSAT, NONUDF 

integers N pb N pc , Np N k> M, Q h Q s , 

Qp> Qn Qf 

real 0 CS , pp p m , dp E, v , G, /, m, w 

real arrays K v( , Kp y d k , ^(1,1000), 

P f (75,1000), P c ( 1,62) 

maximum dimensions E c(t E c p E cm , A cx , C cx , D cx , D*? x , 

A*? x ( 3,3), ap 

a €> N cT t X> M dT ( X' 

^ cM f X ♦ M cM t X > ^ csz ’ L jc (1,6), 

M CX ,N CX (3,7V fc ), £>„(10,6), AINF (6,1000,8), 

( K)p,s> (Qp,s> 

( ty)p,s , (U000) 
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string arrays title (80 characters) read in. 

P$ (eight spaces per field, N p ? fields) 

C e \ (six spaces per field, six fields) 

C € 2 (six spaces per field, six fields) 

P cp (six spaces per field, N pc fields) 
codes C 

current dimensions N p ^ N pc , Np N ms 

real arrays K v %, Kp, dec* P((75,N() 

current dimensions P c (\,N pc ), 

AINF (6 ,N fr 8), \ XtP ,s(UN ( ) t 
?x,y,P,sS l >N() 

The subroutines are as follows: 

INVA inverse of an array 

GACD3 generates composite three-dimensional elastic and thermal properties and 

the two-dimensional thermal properties 
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BLOCK DATA DISP (String) and RESF (String) 

GPCFD2 generates composite two-dimensional elastic constants and constitutive 

equations 

COMSA generates the ply strain and stress states due to applied loads, checks for 

ply failure and interply delamination, and generates the ply stress and 
strain influence coefficients 

INHYB generates ply level properties with the aid of subroutines FIBMT, HTM, 

COMPP, and FLEXX 

BANKRD/IDGER generates constituent properties by using the data base 

FBMTDATA.BANK and arranges them in a proper format so as to input 
to INHYD 

FESTRE computes Poisson’s ratio mismatch between the plies and the composite 

EDGSTR computes interlaminar free edge stresses 

MCRSTR/MINCOF generates the microstresses and the corresponding influence coefficients 
STRCNF generates the stress concentration factors around a circular hole 

NUDIFS generates the Poisson’s ratio differences within the plies and the probable 

locations of delamination around the free edge of a circular hole 

MSCBFL performs failure load analysis based on first ply failure/maximum-stress 

criteria and prints the summary 

AMINF minimum value of an array 

AMAXF maximum value of an array 

FLRLD determines the failure load, failure mode, and the ply location 

These subroutines are described in detail in the next section. 

INPUT title, Nfy N k> N ms , CSANB, BIDE, RINDV, 

COMSAT, NONUDF, t b 9 b 7), M b fiber 
name, matrix name, k vb kj^ k SCJ N cx , M cx , 

DM CX , P u , P ( 

(For substitution and definition, see appendix A.) 

Subroutine Description 

Subroutine INVA (N,A,C). — This subroutine computes the inverse of a square matrix A by Gauss 
elimination and stores it in array C. The check 

|A|*0 

is made and, if satisfied, the program continues; otherwise, the message SINGULAR MATRIX is 
displayed. The subroutine inputs are N, the matrix order, and the matrix A. The output is 

A -1 — C 

Subroutine GACD3(C). — This subroutine generates the three-dimensional hygrothermoelastic 
properties of the composite about its structural ( x,y,z ) and material (1,2,3) axes. The angle 6 is 
measured from x of the structural axes system. (See fig. 6.) In figure 6, replace xy etc. by 1 1 etc. and 
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Figure 6.— Ply orientation geometry. Composite structural axes, x,y,z ; composite material axes, 1,2,3; ply material axes 
(coincides with fiber direction), lj-,2,-,3,-’ 


measure 0 from the material axes to obtain properties about the material axes. These composite 
properties are generated from the following equations: 

r n , n , , -j 

[£J=r E (z e+ i-z (i )[R l ,] T [E (i ][R ei ]+ L HjiSj] 

L /= 1 7 = 1 J 

N t 

Kj = - IE C ) E (««+ 1 

<c (=1 

1 N< 

Wc\= 7 I £ V] E + 

'c /=] 

The arrays \a c l [0 C } 9 and in the preceding equations are given by 

W = [°icxx (x cyy c Xczz Gi cyz c *czx Gi cxy] 
l@c\ “ [&cxx @cyy @czz &cyz $czx$cxy] 
and 


f«fti = [«n l«f22«03 0 0 0 ] r 

(/?&'! = [0fll/3f220f33 0 0 0] 1 * * * * * 7 ' 

For all practical purposes, the two-dimensional thermal coefficients of expansion about the 
composite structural axes are the same as a cxx , a C y y , and a cxv in the array |a c j for the three- 
dimensional case. 
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The matrix [ZTJ 1 is given by 



Note that for the case of an anisotropic material, the elements (1,6), (2,6), (3,6), and (4 5) and their 
symmetric parts will not be zero. 

The matrices [£#] 1 and [/?#] ~ * are given by 



20 



[*<>]= 


cos 2 0 

sin 2 6 

0 

0 

0 

^ sin 26 
2 

sin 2 0 

cos 2 6 

0 

0 

0 

- ]r sin 26 
2 

0 

0 

1 

0 

0 

0 

0 

0 

0 

cos 6 

sin 6 

0 

0 

0 

0 

-sin 6 

cos 0 

0 

■sin 26 

sin 26 

0 

0 

0 

cos 26 


where 0 = 6 (i for properties about the composite material and d = d (i + 6 cs for properties about the 
composite structural axes. (See fig. 6.) 

The matrix [S,y] is given by 


A2 -A2 0 0 


0 -AB 


-A2 A2 0 


0 0 AB 



0 

0 

0 


0 0 
0 0 
0 0 


0 0 0 
0 0 0 
0 0 0 


- AB AB 0 0 


0 B2 I . 


Here A = sin 20, -sin 20, and B = cos 20, -cos 20, where / >1 and denotes the ply index. 

The angles 0, and 0,_j (fig. 6) are given by 

~ ®fi ®C5 
0/-1 + ®cs 

The composite heat capacity is the same for both the two- and the three-dimensional cases. It is given by 

1 N ‘ 

h c - ~ I] h a t(j 

‘c i = 1 


and t c . is given by 

l c~ l e 
i= l 


The composite three-dimensional heat conductivities along the composite material axes, assuming an 
orthotropic composite, are given by 


*di - : £ hi^nx cos2 fy+ A f22 sin 2 Og)j 
l c i- 1 

1 Nt 

K C 22 = 7 L s * n2 K(22 COS 2 0g)/ 

*c /- 1 

a c33 * c /- 1 x A f33// 


The angle is measured from the material axes (fig. 6) 

The composite two-dimensional heat conductivities along the composite structural axes are given 
by (see ref. 9 for the transformation equations) 


K C xx- 7 2^ *&(Km cos2 ^ + A f22 sin 2 0)i 
l c /-I 

1 N * 

K cyy= 7 £ sin2 0 + A f22 cos2 #)/ 

/=1 

1 

A CW = ^CXy = 7" 2^ *ff(^f22“ A m)/‘ s * n 2®| 
tr i — 1 


Kczz-Kcx 


The angle 0 in the last set of equations is measured from the composite structural axes and is equal to 
+ The in PUts to the subroutine are N ( , z &+ ,, z ff , 0„, 0 P/ , [£>], //,, (a ff J, A ffj and {K e }, which are 
all global. The variable N t is input data. The remaining quantities are either generated or are 
transferred from information stored in PL(11,1), PL(13,I), PL(15,I-23,I), PL(8,I), PL(24,I) to 
PL(26,I), PL(30,I), PL(27,I), and PL(29,I). The outputs are f c and the arrays are [ isj -1 , [a c J, [E c \, 
h c , and [K c ). The composite thickness t c is stored in PC(2). The arrays [fj - [a c ), and [fj for both 

composite material and structural axes are printed out under the headings 3-D COMPOSITE 
STRAIN STRESS TEMPERATURE MOISTURE RELATIONS-STRUCTURAL AXES and 3-D 
COMPOSITE STRESS STRAIN RELATIONS-STRUCTURAL AXES. 

The composite material axes properties [E c \ and (a c ] are stored in PC(3) to PC(14) as global 
variables. The corresponding moduli are stored in PC(19) to PC(30). The three-dimensional heat 
conductivities and heat capacity along the material axes are stored in PC(15) to PC(18). The two- 
dimensional thermal coefficients of expansion along the structural axes are stored in PC(48) to 
PC(50). The two-dimensional heat conductivities and heat capacity along the structural axes are 
stored in PC(51) to PC(54). Note that the heat capacity is a scalar quantity and is independent of the 
reference axes. Therefore, PC(54) equals PC(18). The moisture diffusivities and expansion 
coefficients are stored in entries PC(55) to PC(62). 

Subroutine BLOCK DA TA . — In this block, the strings C e \ and C e2 , which are printed out with the 
composite constitutive equations, are defined. The string C el contains the resultant force notation 
N cx> N cy , N cxy , M cx , M cy , and M cxy . The string C e2 contains the notation for the corresponding 
displacements. 

Subroutine GPCFD2 (RESF, DISP, PROPC ). — This subroutine generates the required section 
properties and the force-deformation temperature-moisture relations for a two-dimensional 


multilayered composite. It also generates the plane-stress elastic constants for the composite. The 
force-deformation temperature-moisture relations generated in this procedure are defined in the 
following equation: 


(VcJ 


IA CX ] [C cx ] 


t csx 


l^cr,jr] ( 

1 

l^cM/xl 

{M cx } 


I^cd \P C A 


f w cbx , 

1 i 

< M cT,x\ 

1 1 

[ l^cM ( x\ 


The generic equations for the elements in the arrays [A cx ], [ C cx ], [D cx \, {N cT(X }> \M c j iX j, 
and [M cM(X J are 

N t N f - 1 

IA CX ] = ^ (z$ + \ — Zfi)[R&] T [E$] “ 1 [/?#•] + S El j[Sj\ 

/= 1 j = l 

N t N ( - 1 

[c cx ] = \ E (4+i-«i)W*] r [^]-W+ E z rpj Hj[Sj] 

L i= i y = i 

N, N t - 1 

[£> c J=i E (4+l-4)[^d r I^«] - W+5 E^/W 

J / - 1 Z 7 = 1 

N t 

l N cT t x\ = E ^T a (Zfi + j -2 e )[/? e ][£'&] ~ 1 fotff] 

1 = 1 

N, 

i n cm,x\ = E (za + 1 - - 1 

/■= 1 

1 

[A/ c rJ = , E &T s (zl + \-zl)[R(,\ T [E b ] - 1 {offtj 
z /=! 

1 Nf 

[aw=, E M e (4 +1 -zi)[/? ff ]n£ff]-‘(w 

2 /=l 

where A T & =T & — T CU{ 


The arrays [or ff ], (/3 e j, [/? a ], [£ff], and [Sy] are 


j3 

II 

J3 

! 11 a 22 

0 ]J 


- W[ | &22 

0 \J 



cos 2 0 

sin 2 0 

]r sin 20 
2 

[*ff] = 

sin 2 0 

cos 2 0 

- ^ sin 20 


-sin 20 

sin 20 

cos 20 


23 



0 


J_ v (2i 

E n 1 e (22 


m = 


^ — 0 

E m E ai 


0 


0 


1 

G n2 i 


5/22 -Sj\\- 


1 

4 


(sin 26 j- sin 20/ _ j) 2 


5/21 -5y 12 - ~5yii 

5/32 = 5/23== ^(sin 20/ -sin 20 / _ l )(cos 20/ -cos 20 / _ 1 ) 

5/31 = 5 y 1 3 = — 5y 2 3 

S/ 3 3 = ^(COS 20/ -cos 20 / _ l ) 2 


Here 0/ equals the 0 cs + fy (fig. 6). The reduced bending rigidities (ref. 6) are generated in this 
procedure according to the equation 

D& = [Doc “ C cx A cx E cx\ 

The reduced axial stiffnesses are generated in the procedure according to the equation 


a r — \ a _ c n ~ 1 

^ CX — t^ 1 CX ^cx^cx 


^cxl 


The two-dimensional composite elastic constants are generated from the following equation 
(assuming T (i - 7} for /= 1 to N ( and M & = M? for / = 1 to N$: 


J N t 

[E cx ] “ 1 = t ( E (Z&+1- T [E&\ 

\ / = 1 


N, 


[/?«] + E «/[Syl 

7=1 


> 


where 

A = ^ % 

/=! 


The inputs to this subroutine are M#, 0/ (relative to composite structural axes), Hj , and the 

ply elastic constants. These quantities are global and are located, respectively, in PL(7,I), PL(50,I), 
PL(72,I), PL(14,I), PL(9,I), and PL(31,l)to PL(42,I). The arrays [/y r , [f#]' 1 , [/?*], and [Sj] and 
the dimensions zq are generated within this subroutine. 

The outputs are the force-deformation temperature-moisture relations, which are stored in tne 
global arrays ACX=A cx > RAC-v4^., CP C = C cx , FL X=D cx , RDC =D^ X , NSDT = N c j (X ^ 
MSDT =A/ c7 '^, NSDH = and MSDH=M cA * t x- These are printed out under the heading 
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FORCES FORCE DISPLACEMENT RELATIONS DISPL T-FORCES H-FORCES. The 
reduced bending rigidities are printed out under the heading REDUCED BENDING RIGIDITIES. 
The reduced axial stiffnesses are printed out under the heading REDUCED STIFFNESS MATRIX. 
The inverse of the constitutive equations 


Kvl 

ICcx ] 

[C cx ] 

[D cx ] 


are printed out under the heading DISP DISPLACEMENT FORCE RELATIONS FORCES. 
The distances z c > and z& are stored in PC(31,I), PL(10,I), and PL( 1 1 ,1), respectively. The two- 
dimensional composite stress-strain relations are stored in PC(33) to PC(38), and the two- 
dimensional composite moduli and Poisson’s ratios are stored in PC(39) to PC(47). The two- 
dimensional thermal properties are stored in PC(48) to PC(54), as is described in the section 
subroutine GACD3. 

Subroutine COMSA (M ). — In this subroutine the stress and strain states of each ply are computed 
given the edge membrane forces, the ply temperature, and the changes in curvature. In addition, two- 
ply, combined stress-strength criteria and the interply delamination criterion are generated. Also 
generated are the ply stress-strain influence coefficients. The equations programmed for the it h strain 
and stress states are 


iff,! = [R a ][A cx ] - 1 (\N CX \ + \N cT(X \ + \N cM(X \ +[C cx }\w chx \) 

M = [E*] ~ 1 [R ( ,\[A CX \ ~ > (\N CX \ + \N cTi J + f N cM(X \ + [C cx ][w c/)V ]) 

~ [^f,l _ 1 (^f/!«f/i + +2[/?ft][w' c . /)V |^ 


The reference plane strains e csx and the curvature changes are computed from 


(Cxvi 



[Q-.vl 

- r 



K-a-I 

M 

\ N cT ( x\ | 

♦ 

\N c m (X \ 


1 w cbx^ 

i 

[Qvvl 

[D cx ] 


i 

fi 

1 1 Mcx \ 1 

1 i 

| \M c Tfx\ J 

1 


l| 


when either the membrane force or the moments or both are given. 

The strains are generated locally in EPSL and SIGL, respectively, and are stored in PL(64,I) to 
PL(69,I). The matrices [R (i \ and [£>,•] are generated locally from information transferred from 
PL(14,I) and PL(31,I) to PL(42,I). The distance z & , the ply temperature T &9 and the ply moisture M ti 
are transferred from PL(1 1,1), PL(50,I), and PL(72,I), respectively. The remaining matrices are 


A t v 

— 

ACX 

c cx 

— 

CPC 

N cT,X 

- 

NSDT 

N cM(X 

— 

NSDH 

"ex 

— 

NSB m 

M cT,X 

— 

MSDT 

Mcm t X 

— 

MSDH 

M cx 

— 

MSB,„ 


and w cbx -+ WXX m (local curvature from bending analysis), where m denotes the load condition. 

It is important to note that the stress analysis in the coded form also handles the case where both 
the reference plane membrane strains and the local curvatures are given. In this case the ply strains 
are given by 


25 



l e cxi I ~ csx 1 ^ 1 


where [e cxi ] is the /th ply strain along the structural axis, \e csx \ is the reference plane membrane strain, 
Z is the distance from the reference plane to the centroid of the /th ply, and \w cbx \ is the local 
curvature. These variables are read in the array D nrr where m denotes the load condition. 

The corresponding / th ply stresses are given by 

[cr^j = [£fl] * {[^y] ^cxi ~ ^ y 1 “ 

^ 'Tft = — T cu j 


where [o lT is the /th ply stress along the material axes, [E e j is the /th ply elastic constant about the 
material axes, R#. is the transformation matrix of the /th ply, (e c */j is the /th ply strain along the 
structural axes as given by a previous equation, T & is the temperature of the /th ply, T cui is the cure 
temperature of the ith ply, [a & ] is the thermal coefficient of expansion of the /th ply along the 
material axes, M# is the moisture content of the /th ply, and (£?#) is the moisture expansion coefficient 
of the /th ply along the material axes. 

The displacement force relations are printed out in the following format: 


DISPLACEMENT DISPLACEMENT FORCE RELATIONS FORCES 


I u cx) 


f Wr 


Wed [O c J 

[C cx ] [D cx ] 


Wcx) 


iM cx \ 


Two similar sets are printed out. In the first set, the displacement and force vectors are in symbolic 
form. In the second set, the displacement and force vectors have their numerical values. (See outputs 
of trial cases, app. B.) 

The failure criterion may be determined by either of the following methods: 

(1) Modified distortion energy 


F= 1 - 



S (22 



PL(62,I) 


The parameters a and 0 are specified as follows: 


a = 
0 = 



a n{ >0 
°n \ <0 


°n2 — 0 
o?22<® 


3 nicC 


(min (5 fllc , SniCZ>) 




■(; 


} mr 

'me 


0=T 
0 = C 


a=T 


a = C 
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Kn2a&-Ki\2a0‘ 


(1 +4^12" ^ 13)^22 + U ~ y f23)^ni 
Eg\ l^f22(2 + 12 + ^fl3)(2 + ^f21 + ^23)) 1 /2 


I BET(1, 7) a, (3=T 

BET(2, 7) a = C, 0=T 

BET(1, 8) a=T, 0 = C 

BET(2, 8) a,$ = C 


The multiplier of K?\2al3 was generated in the main program and is stored in PL(61,I). The constant 
K'xiap constitute theory-experiment correlation factors. These are set as unity in COMSA. However, 
the user can modify the correlation factors if he/she wishes, by redefining the matrix BET in the 
subroutine COMSA. 

(2) Hoffman’s criterion (ref. 9) 


Snic- min ( 5 mo Snicz?) 

°}[\ - a n\ a P22 , <$. 22 , , Smc~ Smr „ , ^12 

— e c + c c + “c c °n 1 + e a f22 + -r — 

^nic^niT ^nic^fiir ^220^221 sf ]2 s 

F>0 no failure 
F = 0 incipient failure 
F<0 failure 

The interply delamination criterion for the yth interply layer at the mth load condition is governed by 



>PL(71,I) 


[l - (J^) J —PL(63,I) when />1 

A <pj= ^ ( e C yy “ e cjcx-)(si n 20,- -sin 2 0/_,)+ ~e cxy (cos 26 ; -cos 26 j)(e c J = {A cx \ ~ 1 

K*] "b “i" [QJ ! ^cbx)) 


or by the displacement force equation described previously. 

The inputs to the subroutine are the ply angle measured from the structural axes (0/, from 
PL(14,I)); the distance from the reference plane to the centroid of the ply (z 6 , from PL(1 1,1)); the ply 
temperature (T &t from PL(50,I)); the interply delamination limit from PL(60,I)); the ply 

thermoelastic properties stored in PL(24 to 26,1) and PL(31 to 42,1); the ply extensional and coupling 
rigidities, A cx = ACX and C cx -C PC; the local curvatures w c ^=WXX; the adjustment constants 
*n 27 T=BEf(l, 7), ^ 12 C 7 = BET (2, 7), A^n 2 7 ’c = BET(l, 8), and ^ 12CC =BE T (2, 8); and the load 
conditions N cx = NBS(m). 

The subroutine outputs are the modified distortion energy PL(62,I), Hoffman’s criterion PL(71 ,1), 
the interply delamination criterion PL(63,I), and the adjacent ply relative rotation (A <pj> from 
PL(70,I)). 

Subroutine EDGSTR . — This subroutine computes the interlaminar stresses o zv o zyi and o ^ near a 
straight free edge region of a finite width, infinitely long plate under uniform extension. The 
equations used are based on an approximate formulation analogous to that in reference 18. The 
calculations are performed in two parts. The first part consists of computations of decay lengths for 
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the interlaminar stresses. The decay length is a measure of a free edge region in which the 
interlaminar stresses may be significant. This is achieved in the main program. The second part uses 
this information to compute the interlaminar stresses in the subroutine EDGSTR. The pertinent 
equations are discussed in the following paragraphs. Note that in the case of hybrid composite plies, 
the calculations are repeated not only for the primary composite but also for the secondary composite 
by using the appropriate ply constituent properties. The primary and the secondary composites are 
distinguished by using the letters P and S, respectively, in the Fortran variables. In the case of biaxial 
loading, this subroutine is bypassed as there are no free edges. 

Part 1 . — Decay length or boundary layer width computations . The interlaminar stresses near the 
free edge are assumed to decay exponentially. The decay length is calculated with the aid of the 
following equations: 



where 

« = &' ( 0 . 001 ) 


and 


[X| 


Qn -f * 

Efyy L V 4(1 -k v )kf 


) 


1/2 


The calculations are repeated for each layer. Quantities and X, are stored in arrays YPL and 
PLMDAY. These quantities pertain to the free edge parallel to the load axis X . The corresponding 
quantities for the load axis parallel to Y are stored in arrays XPL and PLMDAX. These are 
computed by replacing E^ v with E ^ in the preceding equations. For the intraply hybrid composite, 
the respective arrays for the secondary composite are denoted by YSL, SLMDAY, XSL, and 
SLMDAX. Note that the letter P is replaced by S. This notation is followed consistently throughout 
the text. The labeled common block ILAB6 is used to store and pass these data to subroutine 
EDGSTR. 

Part 2. — Interlaminar stress computations . In the EDGSTR subroutine, the ply stresses PL(67,I) 
to PL(69,I) are transformed to the structural coordinate system jc, y, and z. These stresses are stored 
in the matrix SIGMA (3,1) for each layer. The interlaminar stresses [a fez J are computed with the aid of 
the following relations: 


/ 



for / = Nf_i to Ay 2 + 1 


No 


2 



2 


Nf 

°fyy 

2 


The interlaminar shear stresses and [o ^ x j are calculated by 


j + 1 1 . . 

^ c^fyyl\ 

i _ a j = Nt 

°&y - 1 ) 

ft 


fori = N ( to y + 1 
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and 



for i = NAo 




+ 1 


In these equations, the computations are started from the top layer (/ = N f ). After the midplane is 
approached (i = N/ 2+ 1), the calculations are repeated starting from the bottom layer (/= 1) and 
continued until i becomes (N/ - 1). 

The interlaminar stresses are stored in the arrays YSZZP, SZYP, and SZXP for the primary 
composite and in the arrays YSZZS, SZYS, and SZXS for the secondary composite. They are, 
however, made dimensionless by dividing by the applied normal stress o txx . 

Subroutine STRCNF . — This subroutine calculates the stress concentration factors around a 
circular hole due to membrane loading. The equations used are taken from reference 19 and are 
strictly applicable for infinite plates. Three factors are computed in the subroutine and are defined by 
the following equations: 


K 


I XV 


a ee 

°X\'Oo 


K lyy = 


°ee 

°yy o° 


K 


\xy 


a 66 

°xyoo 


Quantities a vvoo , a yyoo9 and a xyoo are the applied stresses, and is the hoop stress at any angle 0 from 
the load axis. The stress concentration factors are stored in the local arrays XK1, XK3, and TEMP. 
The expressions for K Xxx , A' lvv , and K\ xy are the following: 



In the preceding expressions, Ecu and E crr are the composite moduli in the tangential and radial 
directions at angle 0. Angle 6 is measured from the x-axis for K Xxx and K Xxy and from the>>-axis for 
K ]yy . The program rearranges the computed K Xyy values so that they correspond to the same location 
as those of K Xxx and K Xxv . 

Subroutine NUDIFS . — In this subroutine, the Poisson's ratio differences between the adjacent 
plies and the composite are computed around a circular hole at 5° intervals. The products of the 
differences and the corresponding stress concentration factors are computed next. These products are 
expected to provide insight into the probable delamination locations. It is assumed that onset of 



delamination is likely to occur at the locations for which the product of Poisson’s ratio mismatch 
with the corresponding stress concentration factor is a maximum. Accordingly, these products are 
computed at 5° intervals and the maxima are calculated. Two sets of tables are the output from this 
subroutine. The first table comes out optionally if the boolean NONUDF is set to FALSE. It contains 
all the details of the computations. The second table consists of the summary of results, with notes on 
the maxima and the locations. The following are the programmed equations: 

At any angle 9 the Poisson’s ratio is computed by 


v crt err 


^CX2_ _ f 1 + 2v cxy + V os 2 e sin 2 9 

- E CXX ' E CXX E C yy ^CXV ' 


The ply Poisson’s ratio is given by 



1 

E ai 



cos 2 0 sin 2 0 




The difference in Poisson’s ratio between the /th and (/ + l)th plies is given by (v 1 ^ 1 - vjr/), and the 
difference with respect to the composite is given by (u^ { - v crt ). These are stored in the arrays A2 and 
A3, respectively. The products of K ]xx , K ]yy , and K Xxy with A3 are computed next and are stored in 
the arrays A5, A6, and A7, respectively. The maxima and their location in each of the four quadrants 
(0-90, 90-180, 180-270, and 270-0) are computed by calling the subroutine AMAXF for the three 
arrays A5, A6, and A7. The values of stress concentration factors are passed through the labeled 
common block ILAB8 from the subroutine STRCNF. 

Subroutine MSCBFL (AIIWF).— A complete laminate failure stress analysis, based on first-ply 
failure and the maximum strength criteria, is performed in this subroutine. The inputs to this routine 
are the ply allowables S mc , S mr , S G2 o Sf22T> and S n2 s anci the ply influence coefficient matrix 
A1NF. The ply stress allowables are generated by the INHYD routines and are stored in the ply 
properties array PL. These are accessed through the labeled common block ILAB2. The ply stress 
influence coefficients are generated by COMSA and the main program and are passed to the present 
routine by the subroutine argument. 

The failure stress for a particular ply due to a specific loading is given by the ratio of the allowable 
strength to the ply stress influence coefficient. For example, the failure stress due to a tensile load is 
given by 


«_ jfojr 

c Factl' 


where Fact!' is the stress influence coefficient for /th ply due to unit tensile loading, S l n]T j s the 
strength allowable for /th ply in longitudinal tension, and S', is the failure stress for the ith ply due to a 
tensile loading. The failure stresses are stored in the matrix FAILD. In the case of 
temperature/moisture presence, the allowable strengths are updated to take into account temperature 
or moisture stresses; the failure stresses are computed with and without the effects of temperature- 
and moisture-induced stresses for comparison. The program considers primarily five different 
loadings, longitudinal compression and tension, transverse compression, and tension and inplane 
shear. 

After the failure load computations for each ply are determined, the active failure mode and the 
corresponding failure strength for each type of loading are determined by calling the subroutine 
AMINF. This subroutine returns the value of the minimum failure load, the ply number, and the 
failure mode as output. The output from this subroutine is printed under the heading LAMINATE 
FAILURE STRESS ANALYSIS. 

Subroutine MCRSTR .— This subroutine generates the microstresses in the ply constituents due to 
the inplane loading. These are stored in the ply microproperty arrays PLMP and PLMS for the 


primary and the secondary composites. The ply constituent properties and the applied loads are 
inputs to this subroutine. They are accessed with the aid of the common blocks PBANK, MFBANK, 
ILAB2, ILAB5, and ILAB9. The PLMP and PLMS each contain 41 entries which are explained in 
the following list: 


Code 

Algebraic notation 

Fortran variable 

name 



PLM(1,I) 

°m\ \L 

SMIL 

PLM(2,I) 

°m nr 

SM1T 

PLM(3,I) 

°f\\L 

SF1L 

PLM(4,I) 

°f\\T 

SF1T 

PLM(5,I) 

(A) 

22L 

SM2AL 

PLM(6,I) 

JA) 
a m 22 r 

SM2AT 

PLM(7,1) 

JB) 
a m 22L 

SM2BL 

PLM(8,1) 

(B) 

°m 22 r 

SM2BT 

PLM(9,I) 

JB) 

a /22L 

SF2BL 

PLM(10,I) 

(B) 

°/22T 

SF2BT 

PLM(11,I) 

JA) 
°m 33Z. 

SM3AL 

PLM(12,I) 

(A) 

a m 33 r 

SM3AT 

PLM(13,I) 

JB) 
a m 33£ 

SM3BL 

PLM(14,I) 

JB) 
a m 33 T 

SM3BT 

PLM(15,I) 

JB) 
°f 33L 

SF3BL 

PLM(16,1) 

JB) 
°f 33 r 

SF3BT 

PLM(17,1) 

(A) 
a m 12 

SM12A 

PLM(18,I) 

JB) 
°m 12 

SM12B 

PLM(19,I) 

JB) 

a /l2 

SF12B 

PLM(20,I) 

JA) 
°m 13 

SM13A 

PLM(21,I) 

JB) 
a m 13 

SM13B 

PLM(22,I) 

JB) 

ff /13 

SF13B 

PLM(23,I) 

JA) 
a m 23 

SM23A 

PLM(24,I) 

JB) 
a m 23 

SM23B 
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SF23B 


PLM(25,I) 

°f 23 


PLM(26,I) 

°m\\ 


PLM(27,I) 

a f 11 


PLM(28,I) 

JA) j 
°m22 I 


PLM(29,I) 

(B) 1 

a m22 1 

(B) / 

°/22 1 

Microstresses due 
> to temperature 

PLM(30,I) 

gradient A Tt 

PLM(31,I) 

(A) I 

I 


PLM(32,1) 

(B) 
a m 33 


PLM(33,I) 

(B) 

a /33 j 


PLM(34,I) 

1 

a i nil 


PLM(35,I) 

a f 11 


PLM(36,1) 

(A) J 

°m 22 j 


PLM(37,1) 

JB) f 

a m22 \ 

Microstresses due 

PLM(38,I) 

(B) / 

a /22 \ 

to moisture M f 

PLM(39,I) 

(A) 1 

j 


PLM(40,I) 

(B) 
°m 33 


PLM(41,i) 

(B) 

°/33 i 



SM11DT 

SF11DT 

SM2ADT 

SM2BDT 

SF2BDT 

SM3ADT 

SM3BDT 

SF3BDT 

SM11DM 

SF11DM 

SM2ADM 

SM2BDM 

SF2BDM 

SM3ADM 

SM3BDM 

SF3BDM 


In this list, entries 26 to 41 are suppressed automatically if the temperature gradients and moisture 
contents are not present. The superscripts A and B refer to two regions as described in figure 4. 

The microstresses are calculated with the aid of the following equations: (For notation and sign 
conventions, see figs. 4 and 6.) 

Ply microstresses due to a longitudinal stress j are given by 


a m 1 1 — i )<*n 1 


°/\ \ -(£/ii/£Vii) a n! 


(A) , 


, n2 )(E m /E n )a n] 


(B) (B) _ 1 - Vkf 

°m22 ~ a J22 ~ ^J=~ a n I 


(A) _ (A) 
°m 33 “ °m22 



(B) (B) 

CT /w33 — °m22 


(B) (B) 

°m = °fl 2 


Ply microstresses due to a transverse stress a m are given by 

a "'" = V"'" ~E^) m 


_( Em\„ 

< 7/3 1 “ yf\ 2 ~ "02 J a (22 


°m22 = (E m / E 2 )0(22 


0< m22 ~ (Eni/ E 2 ) a n2 


< 7/22 — (E f 22 ^ E 2 )o ^22 


where £2 > s g' ven by 


E 2 = (l-y/k f )E m + 


\lk f E m 




o,nh = ( v m /v m)(E m /En.2( a a2 


(B) 

a /«33 = 


1 -\/k 




^T f 


a f22 


(B) (B) 

ff /33 — < 7/7733 


Ply microstresses due to inplane shear stress (cr n2 ) are given by 


< 7 ml 2 = (G„/G| 2 ) ff fl 2 


< 7 ml 2 = (Gfi 2 / G l 2 )ff(i 2 


< 7/12 = (Gf| 2 / Gl 2 )< 7 n 2 


where G \ 2 is given by 
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(l-^) 


VF/G. 


'-M'-t) 


(A) _ (A) 
°m \3 -°m \2 


(B) _ J B > 

1 3 a ;?7 1 2 


(B) _ (B) 
°yi 3 -°f \2 


Ply microstresses due to through-the-thickness shear stress oq 3 are given by 


a m23 = (G m / Ge23) a m 


a rn23 = ^ G 23 /G m) a m 


where G23 is given by 


M'-fe) 


(B) (B) 

°/23 = CT m23 


Ply microstresses due to temperature gradient A 7 ) are given by 


°m\ 1 -(<*ni ~ a m)^T ( E n 


°fi 1 = («n 1 _ «/i 1 )AT ( Ef\ 1 


°n,22 = (°‘a2- a m) AT l E n 


JB) _J B > 
a m 22 -°f 22 


1 -sTkf 


(A) (A) 

a m33 — a m22 


JB) _J B > 
ff m33 - a m22 
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(B) (B) 

ct /33 = a J22 

AT ( = T ( - T cu 

Ply microstresses due to moisture M f are given by 
<T ml\ = (@n l ~ & m )M ( E m 

a f\ \ =^fi \M ( Ef ] | 

a< m22 = (002 ~ 

(B) _ 1 -V£7 (A) 

a m22 ~ 7=^ a m22 

^k f 

(A) _ (A) 
a mi 3 - °m22 

(B) _ (B) 

°m 33 °m22 


(A) _ (A) 
°f 33 — °>33 


Subroutine MINCOF . — This subroutine generates the microstress influence coefficients for each 
different material system used in the layup. The equations used are similar to those programmed for 
MCRSTR. However, the influence coefficients are based on the application of unit load in a specific 
direction or unit temperature difference or unit moisture content. The influence coefficients are 
stored in the matrix P1NF. This matrix has 17 entries. They are described in the following list: 


Entry 

Algebraic notation 

Fortran variable 

PINF(1 ,K,NLD) 

a m ll 

SMI 1 

PINF(2,K,NLD) 

(A) 

a m22 

SM22A 

PINF(3,K,NLD) 

JB) 
a n 722 

SM22B 

PINF(4,K,NLD) 

(A) 

°m\2 

SM12A 

PINF(5,K,NLD) 

JB) 
a m 12 

SM12B 

PINF(6,K,NLD) 

(A) 
a m 13 

SM13A 

P1NF(7,K,NLD) 

JB) 
a m 13 

SM13B 

PINF(8,K,NLD) 

JA) 
a m 23 

SM23A 

PINF(9,K,NLD) 

JB) 
°m 23 

SM23B 
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PINF(10,K,NLD) 

„< A > 
°m 33 

SM33A 

PINF(1 1 ,K,NLD) 

(B) 
°m 33 

SM33B 

PINF(12,K,NLD) 

°f 11 

SFI 1 

PINF(13,K,NLD) 

JB) 
a f 22 

SF22B 

PINF(14,K,NLD) 

JB) 

<7/33 

SF33B 

PINF(1 5,K,NLD) 

a fl2 

SF12 

PINF(16,K,NLD) 

ct /13 

SFI 3 

PINF(17,K,NLD) 

JB) 

a /23 

SF23B 


The dimension K varies from 1 to NMS, where NMS is the number of material systems. NLD 
varies from 1 to 7. The expression NLD= 1 to 5 refers to unit applied stresses in 11, 22, 12, 13, and 
23, respectively. The expression NLD -6 corresponds to unit temperature loading, and the 
expression NLD -1 corresponds to unit moisture loading. 

The microstress influence coefficients are computed for secondary composites and optionally 
computed for intraply hybrid composites. These are stored in the matrix SINF. 

Subroutines AMAXF , AMINF , LOGO , and L0G02 . — These subroutines perform several 
auxiliary duties. AMAXF finds the maximum value of a one-dimensional array and its location. 
AMINF finds the minimum value of a one-dimensional array and its location. These two subroutines 
are utilized by MSCBFL and NUDIFS in conjunction with searching for failure loads and the 
probable locations of delamination. LOGO is a subroutine to generate the ICAN emblem for the 
output. The description of the material and the structural coordinate system by appropriate figures is 
generated by the subroutine L0G02. 

Subroutine INHYD . — This subroutine generates the composite ply properties, necessary for the 
laminate response analysis. The inputs to this routine are the constituent properties which are 
supplied in the appropriate format by the subroutines IDGED and BANKRD. INHYD calls the 
subroutines FIBMT, COMPP, and HTM to perform the micromechanics analysis, including the 
analysis of hygrothermal effects. The ply properties are stored in the array PROPS, which is accessed 
by the main program through the labeled common block PBANK. INHYD is called once for each 
different material system by the main program. The outputs of INHYD show the properties of the 
fiber, matrix, and composite. 

The fiber and matrix properties for the primary composite are read in from the input provided by 
IDGER. These are stored in arrays PF and PM. Similarly, the arrays SF and SM are used to store the 
properties of secondary composite constituents if the composite is of the hybrid type. The program 
then checks for temperature and moisture. The properties of the matrix are updated for the presence 
of temperature and moisture. The following are the equations programmed to account for the 
hygrothermal property degradation: 

The wet glass transition temperature is computed from 

T gwr = (0.005M}-.\M t + \)T gdr 

where T gwr is the wet glass transition temperature, T gdr is the dry glass transition temperature for the 
resin matrix, and A/^is the percentage of moisture by weight. 

The reduction factors X mp and X tp are computed from 

X mp ~ v ( T gwr - T u )/{ T g( j r — T q ) 

— l ! X mp 

where T a is the reference temperature (70 °F), and T u is the use temperature. 



The moduli and strengths of the matrix are multiplied by X mp to obtain the new properties for the 
matrix. The density is given by 

Pmw = Pm 3p m k m M/\00 

The thermal properties, such as heat capacity, thermal expansion coefficient, and thermal 
conductivity are multiplied by the second factor X tp to account for the hygrothermal conditioning. 

After the property arrays PF, PM, SF, and SM are properly filled, the program chooses either 
FIBMT or HTM subroutines to perform micromechanics. The subroutine HTM is chosen if 
temperature/moisture effects are to be taken into consideration. Otherwise, FIBMT is chosen for dry 
room temperature property computations. The outputs from these routines are primary and 
secondary composite ply properties. They are stored in the arrays P and S, respectively. These 
properties are made common to subroutine COMPP through the common blocks ILAB1 and 
ILAB3. The subroutine COMPP is called by INHYD for hybrid composites to compute the hybrid 
composite ply properties. These properties are stored in the array H. One of the arrays P, S, or H are 
passed to ICAN via common block PBANK and the array PROPS. For example, if the ply is made of 
100 percent primary composite only, the array PROPS is assigned to have the same entries as P, etc. 

The subroutine INHYD also calls FLEXX, which performs a flexural strength analysis. However, 
these are only for additional information and are not used by ICAN at the present time. 

Subroutine FIBMT (C 9 F, M f VF, VM t VP , KV, IFLAG).— This subroutine generates properties 
of a ply by using the constituent properties which are supplied from the subroutine INHYD. The 
constituent properties are stored in the arrays F and M; F contains the fiber properties, and M 
contains the matrix properties. The composite properties are stored in the array C, which is returned 
to INHYD. The theory behind the programmed equations is discussed in reference 13. The following 
is a description of each entry in the arrays C, F, and M, with the corresponding algebraic notation: 

Composite Properties Array C(I) 


Entry Description Notation 

C(l) elastic moduli An 

C(2) elastic moduli A 22 

C(3) elastic moduli £53 

C(4) shear moduli Cn2 

C(5) shear moduli 

C(6) shear moduli Ai3 

C(7) Poisson’s ratio 

C(8) Poisson’s Ratio %3 

C(9) Poisson’s Ratio ^13 

C(10) thermal expansion coefficient «ni 

C(ll) thermal expansion coefficient a ? 22 

C(12) thermal expansion coefficient 0*33 

C(13) density Pa 

C(14) heat capacity Q 

C(15) heat conductivity An 

C(16) heat conductivity A 22 

C(17) heat conductivity A33 

C( 1 8) strength Alir 

C(19) strength A 1C 

C(20) strength $n2T 

C(21) strength A 22 C 

C(22) strength S n2 

C(23) moisture diffusivity All 

C(24) moisture diffusivity A’ 22 

C(25) moisture diffusivity A33 

C(26) moisture expansion coefficient All 

C(27) moisture expansion coefficient 



C(28) 

moisture expansion coefficient 


C(29) 

flexural moduli 

E n i 

C(30) 

flexural moduli 

E(22 

C(31) 

strengths (flexural) 

s 123 

C(32) 

strengths (flexural) 

S|'l l F 

C(33) 

strengths (flexural) 

S(22F 

C(34) 

strengths (flexural) 

S(12 

C(35) 

ply thickness 

t( 

C(36) 

interply thickness 

6 ( 

C(37) 

interfiber spacing 

Fiber Properties Array 

?>s 

Entry 

Description 

Notation 

F(l) 

elastic moduli 

Ef ii 

F(2) 

elastic moduli 

Ef 22 

F(3) 

shear moduli 

Gf 12 

F(4) 

shear moduli 

Gf 22 

F(5) 

Poisson’s ratio 

*712 

F(6) 

Poisson’s ratio 

"/23 

F(7) 

thermal expansion coefficient 

«/ 11 

F(8) 

thermal expansion coefficient 

“/ 22 

F(9) 

density 

Pf 

F(10) 

number of fibers per end 

N f 

F(1 1) 

fiber diameter 

d f 

F(12) 

heat capacity 

c f 

F(13) 

heat conductivity 


F(14) 

heat conductivity 

Ef 22 

F(15) 

heat conductivity 

E/33 

F(16) 

strength 

SfT 

F(17) 

strength 

Matrix Properties Array 

S/c 

Entry 

Description 

Notation 

M(l) 

elastic modulus 

E m 

M(2) 

shear modulus 

G m 

M(3) 

Poisson’s ratio 

v m 

M(4) 

thermal expansion coefficient 


M(5) 

density 

Pm 

M(6) 

heat capacity 

C m 

M(7) 

heat conductivity 

Km 

M(8) 

strength 

SmT 

M(9) 

strength 

$mC 

M(10) 

strength 

$mS 

M(1 1) 

moisture coefficient 

0m 

M(12) 

diffusivity 



The following are the programmed equations for the entries in array C: 


Normal moduli: 


Ef i ] — k/Ef\ \ + k m E m 

£ _ E m 

G2 \-Jk}(\-E m /E m ) 

Shear moduli: 



Poisson’s Ratio: 

v m = v m + kf( v f\2- v m) 


v tl3 = v n2 



Coefficients of thermal expansion: 

«/i 1 + ^ m [( a A /£ /i 1) “ a f \ il 

«fii = 
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a C22 = a m ( 1 - -Skf) 


+ a m k f 



“33 = a l72 

Density: 

Pt = Pfk/+ p m k m 


Heat capacity: 



Heat conductivities: 

l k/Kf \ i + k m K m 

y/kf 

K n2 = {\-^k f )K m + K m 

'-M'-fe) 

K m= K m 

In the preceding equations, K m should be replaced by 

K n{ Jk v 

*„,-(• ~'Jk v )K m + 

if there are voids. The quantity K v is the void conductivity. 

Strengths: 

s n l T= Sjjikj--^ k m E m /Efx j) 

The longitudinal compressive strength is computed based on three different criteria, rule of mixtures, 
fiber microbuckling, and delamination. The minimum of the three estimates is returned as Sf UC . The 
equations for the three cases are 



5.- nr (rule of mixtures) = Sj c (kf+ k m E m /Ef\ i) 
S uc (delamination) = ( 1 3S n 2 + $mc) 

G m 

Saic (fiber microbuckling) 



The transverse strengths are calculated from 
S a2T = S mT (FACJ/DENOM) 

Sf22C = S mc /DENOM 

^ _ [(/^i ~ 1 + G m / G fi^FzG 2^/775! 

^n2~ q f 

where F, and F 2 are defined by the equations: 



The variable DENOM is a Fortran variable given by 


DENOM = 
where <p is given by 



F i-l 


Vi 


+ 1) + y ('P 


1)2 


The Fortran variable FACT takes the value k m if 1FFAG is unity. Otherwise FACT takes the value 
unity. This variable is introduced to correlate the strengths of HMS and Kevlar fiber composites with 
the experimentally observed values. The main program 1NHYD checks for these fibers and assigns 
the appropriate values for 1FLAG. 1FFAG is set at zero for other fibers. 
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Moisture diffusivities: 


D (\ 1 = k m D m 

■ D f22 = (l 

D (33 = D B2 

Moisture expansion coefficients: 

0(\ 1 = 0m^m E m^ E (] 1 

0(22 = 0m(\ -'Tkj)( 1 +kjV m E^ x / E nx ) 

0(33 = 0(22 

Flexural moduli E nlF ): 

E n\ f = E(\\ 

E Q2F = E d2 


Flexural Strengths: 


0(23 F 


( Fl ~ l+ fy^) F2GmSms 

G m F\ 


$ei2 f- 

Ply thickness: A default value of 0.005 is set for t p. This is overridden by the user specified value 
the I CAN main program. 

Interpiy thickness and interfiber spacing: 





Subroutine HTM (C 9 b\ M , VF, VM, VV, IFLAG). — This subroutine generates the 
hygrothermomechanicai properties based on the theory proposed in reference 15. The subroutine is 
called only if nontrivial entries for the use temperature and the moisture content (T u ^70 F or 
nonzero moisture content) are present. The equations programmed are mostly those discussed in the 
subroutine FIBMT description. Therefore, only the equations which are different are mentioned 
here. 

Moisture expansion coefficients: 


@172-0 ~Skf)0 m 




s Jk/E(22 0 — 'fkj) E m 




Pm - P a.2 


Strengths: 


So.it- 



FACT 


>f22C 


-fe) 




■FACT 


1 - ( 'ikjEpji / C/22) 


s «=( s €)°" 2 


{^jG 

f 12 / G/12) 


-FACT 


w -(c?) 


1 - {ykjGm / Gj2j) 


Sfi2C= 1.5S n 2 


In the preceding equations, FACT is a Fortran variable which is given by 
FACT = b s /bf 

for Kevlar and HMS fibers. For all other fibers FACT= 1. 

Subroutine FLEXX (C ).— The entries C(32) and C(33) of the ply property array C are generated in 
this subroutine. They are, respectively, the longitudinal flexural strength and the transverse flexural 
strength. The longitudinal flexural strength is given by 
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2-5S n i j 



The transverse flexural strength is given by 


2.5Sq2T 



Subroutine COMPP (IPFLAG, ISFLAG ).- This subroutine is called by INHYD to generate the 
properties of a hybrid ply. The equations are based on the theory proposed in reference 13. The 
properties are stored in the array H. The entries are, however, the same as those of array C given in 
the description for subroutine FIBMT. The inputs to this routine are the primary composite 
properties array P, the secondary composites property array S, and the percentage of the secondary 
composite k sc . The equations are the following: 


Elastic normal moduli: 

£ n ,(H) = E n ,(P) + [E n ,(S) - E n ,(P)]Ar, c 

g t22(P) 

2?r22(H) — 1 + A: ic .[£’f22(P)/ £(22(S) — 1] 
£>33(H) = E m ( P) + [£,33(S) - £ f3 3(P)Rc 


Shear moduli: 


G f23 (H) = 


3 (P) 



Gt?3(PA 

Cf23(S)/ 


Cq 2 (P) 

G n2 (H) 

^_G n2 (P)\ 

1 M 1 G n2 ( ?)) 

G n3 (H) = C n3 (P) + MGn 3 (S)-G fl3 (P)] 
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Poisson’s ratios: 


^ fl2(H) — ^12(P) + ^ 5 c[^fl 2 (S) ~~ v n 2(P)] 


^fl3(H) = v^\i{P) 4- 


* S r[^2(P)-*'fl 2 (S)] 


(l-Ar 5C .)[£ B3 (P)/£ n3 (S)]-Ar 5l 

‘'(’23(H) = v m^) + ^sc-["f23(S) - «'f23(S)] 

Coefficients of thermal expansion: 


«fii(H) = ■ 


«(33(H) = 


«ni(P) + ^5cl[ a fl l(S)£f. l(S)/ E n 1 (P)] an i(P)l 

£ - i» n3 (H)£ 0 3(H)o fl i(H) + (1 - k sc )E 03 (P) [ 

) 


£ 03 (H) 

+ k sc E m ( S) j^af22(S) + ‘ , n 3 (S)af) j(S) 


a ( 22 (H) = (1 -k sc ) 


4" kr 


^a C2 (S) 1 +«'f23( p )J + "PI2( p ) Q! fll] 

j“f22(S) +*'f2 3 (S)j + J'fi2(S)apii(S)^j 


— efi2(H)af] |(H) — ef2 3 (H)a(3 3 (H) 


Density: 

Pf(H) = (l — k sc )pf(P) + k sc p((S) 

Heat capacity: 


«( 2 2 ( p ) + l 'fi 3 ( p ) a ni( p ) 


Cf{H) = jo -k sc )\Cj(P)kj(P)pj{P) + C m (P)k m (P) Pm (P)\ 
+ k sc \Cj<S)kj(S)pj{S) + C m (S)k m (S)p m (S)\ 

+ \k/,P)k„(P) + k sc k v (S)\Mp mS fC ms i^ / pf(H) 
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where p ms{ and C mst are the moisture density and heat capacity, respectively. 


Heat conductivities: 


K n ,(H) = (1 - **)|*/<P)*fli(P) + k m (V)K m (P)) + k sc \k/S)K„ ,(S) + k m (S)K m (S)} 

„ , D , _ (1 - sfWWmP) + ^¥)K m (P) 

A n2^) - , . — , 

1 - V*y<P)[l -K m {?)/K n2 (V)] 

* (S) = ( 1 ~ ^^,(S) + vT 7S)A„,(S) 
t22 1 - VtyS)[l - K m {S)/K n2 (S)} 

»)= ^ ?(P) 

\-k sc [\-K m (S)/K m i?)) 


K m (H) = K a2 (H) 


The void conductivity K„ with moisture content M is given by A:„ = MA r mst . If there are voids in the 
primary composite, A^fP) in the preceding equations for heat conductivities is replaced by 


/f m (P) = [i-V^xF)]/f m (P) + 


v'W m (P) 

[1 -V£7P)][1 -Af m (P)//g 


Similarly, for the secondary composite, K m ( S) is replaced by 


/f OT (S) = [l-V^(S)]A: OT (S) + 


[1 - V^xsj][i 


Strengths .— The longitudinal strengths are based on the rule of mixtures: 


S n]T {H) = S ni y<P)( 1 - k sc ) + S n i j{S)k sc 
S ( i ic(H) = Sn ic( p )( 1 - ^ic) + s n ic(S)A: 5C 

The following are a few intermediate variables defined for convenience in the evaluation of 
transverse strengths: 

Q p = 1 - 2VA: ( ,(P)/7r [ 1 - 2V*/P)/t] 


£> s =l - 2VA-„(S)/7r fl — 2v/A-/S)/ir] 



S mC =rnin[5 mC (P) and S mC (S)] 

S m7 -=min[S mr (P) and S m7 {S)] 

S„, s = rnin[S mS (P) and 5 mS (S)] 

FACT 1 = k m ( P) (for HMS and Kevlar fibers) 

FACT 2 = k m (S) (for HMS and Kevlar fibers) 

FACT 1 = FACT 2 = 1 (for all other fibers) 
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The transverse and the shear strengths of hybrid composites are given by 

/U J (1 -* iC )FACTl A\.,FACT2 1 

5q2t<H)-£C 2(H) |_ £’ f22 (P) D ENOMP + £' f22 (S)DENOMS J Sm7 ' 


Sa2(i^)-E!22{^) L£ Q2 (P)DENOMP + £f22(S)DENOMS J ^ mC 


5 n2 ( H) = 


(1 -kM 
2G nl \ G «2(P) 


G n2 (S) 


FACT2 


['-V¥(-|j§)]) 


-VtyS)(i 


M s) 

G/i 2 (S)> 


Moisture diffusivity: 


D'm= 


D'A S) = 


l-V^XP) + Ar y (P) 
[l-V^fPj^P) 

1 - V^XS) + £„(S) 


£> ni (H) = (1 - k sc )k m (P)Dtf>) + k sc Jc m (S)D’i(S) 


£>e 2 (H) = (l -* JC )[l -2VfyP)]z^P) + *Jl -V^<S)]o^S) 


D m (H)=D a2 (H) 


Moisture expansion coefficients: 


A: m (P)g m (P)£ w (P) 
k f (P)E fn (P) + k m (P)E m (P) 

k m (S)& m <$)E m (S) 
k f (S)E /u (S) + k m (S)E m (S) 


kAP)k m (P)E ni (P) 


^/ii(P) + Ar OT (Pfc(P)-£ /11 (P)] 
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ta®-<us>[' ■ -m] (> + ^, (S) *tXrs!-% l (s)l ) 


0ni(H) = 


[(1 - k sc )k m (P)(3 m (P)E m (P) + k sc k m (S)(3 m (S)E m (S)] 

£ m < H ) 


( 8 ( 33 (H) = ^— j' l >i 3 (H)£'( 33 (H)( 3 fi ](H) + (1 — Ar sc )£'f33(P)[/3(22(P) + |, ?13(P)Gfl l(P)J + ^.rc^t33(S)[£G2(S) 
+ , n 3 (S)/3 m (S)|]A,3 3 (H) 

^(P) = "m(P) + ^/P)t"/12(P) - "m(P)l 

"fl2(S) = »V,(S) + fy{S)[j>/i2(S) - »'/n(S)| 

^P) = "m(P) + M P )I ( ’/23(P) - "m(P)l 
«'o2(S) = MS) + */s)[^ 3 (S) - ^(S)l 

( 3 f22 (H) = (l - Ar sc )[/3f2 2 (P)(l +J' 62 (P)) + » , fl2(P)i8fll(P)J +A: sc[^f22(S)(l +»'<32(S)) 

+ ^l2(SWfi 1(S)| - i' fl2 (H)/3fli(H) - ^(23(H)/3f33(H) 

Flexural moduli: 


£ fll/ KH) = £ m (H) 

£ E2f <H) = £ Q2 (H) 

Flexural strengths: 


Id -MQ 


SC^P 


^f23f( H )- 


2G m (H)} G Q3(P) 


f ^scQs 

G a 3 (S) 


[-V¥0 


-M'-m 


V 

*m(\ 

7T \ 

(1- 

Gm(P) 

V 

G i23 (P) 

G m ( S) \1 


Gm( P) 
c C3 (P) 


)] 


Ga 3 (S)/J 


SJS) 


Sn 25 s(H) = 1.5S fl2 (H) 



Fiber volume ratio: 


kj(H) = */P) + k sc \k/S) - kj{ P)J 


Subroutines BANKRD and IDGER . — These two subroutines do preprocessing to generate 
compatible input data to the subroutine INHYD. The subroutine BANKRD is called first by the 
ICAN main program. The input to this routine is primarily the data supplied on the material card 
MATCRD by the user. These cards indicate the coded names for the fiber and matrix, the volume 
ratios of primary and secondary composites, and their respective fiber, and the matrix and void 
volume ratios. The subroutine BANKRD has its own data base containing the properties of fibers 
and matrices of commonly used materials. This data base is assigned to input unit 8. It is named 
FBMTDATA.BANK. The output of BANKRD are the arrays PFP, PFS, PMP, and PMS. The 
entries in PFP and PFS are the fiber properties of primary and secondary composites. The entries in 
PMP and PMS are the matrix properties of primary and secondary composites. These arrays are 
made common to the main program and the subroutine IDGER through the labeled common block 
MFBANK. The entries of PF and PM arrays are explained in the following list: 


Fiber Property Arrays PFP and PFS 


Entry 

Description 

Notation 

1 

not used 

.... 

2 

fiber density 

Pf 

3 

normal moduli 

Ef\\ 

4 

normal moduli 

Efi 2 

5 

Poisson’s ratio 

v f\2 

6 

Poisson’s ratio 

v f2l 

7 

shear moduli 

Gf\2 

8 

shear moduli 

Of 23 

9 

thermal expansion coefficient 

0/11 

10 

thermal expansion coefficient 

«/22 

11 

heat conductivity 

Kfll 

12 

heat conductivity 

Kf 22 

13 

heat capacity 

C f 

14 

strengths 

Sfr 

15 

strengths 

S fC 

16 

not used 


17 

not used 



18 

not used 



19 

not used 



20 

number of fibers per end 

N f 

21 

fiber diameters 

4 
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Matrix Property Arrays PMP and PMS 


Entry Description 


Notation 


1 

not used 


2 

density 

Pm 

3 

normal modulus 

E m 

4 

Poisson’s ratio 

v m 

5 

coefficient of thermal expansion 

a m 

6 

heat conductivity 

K,n 

7 

heat capacity 

c m 

8 

tensile strength 

S mT 

9 

compressive strength 

S mC 

10 

shear strength 

SmS 

11 

allowable tensile strain 

e mT 

12 

allowable compressive strain 

t mC 

13 

allowable shear strain 

e mS 

14 

allowable torsional strain 

e mTOR 

15 

void conductivity 


16 

glass transition temperature 

Tgdr 


The coded names for the fiber and matrix are stored in the matrix CODES by the main program. 
The entries in CODES are explained as follows: 


CODES(l, 1,1) 
CODES(l,2,I) 
CODES(2,l,I) 
CODES(2,2,I) 


coded name of primary fiber 
coded name of primary matrix 
coded name of secondary fiber 
coded name of secondary matrix 


The subroutine IDGER takes the information generated by BANKRD and arranges it in a proper 
format for the subroutine 1NHYD. These data are transferred to input unit 7 prior to calling 
INHYD. These data are purged at the end of the program execution. 


Data Base FBMTDATA.BANK. 

The constituent properties data base is a unique feature of the computer code ICAN. Its primary 
aim is to reduce the burden on the user by preparing properly formatted data for the program. The 
user only needs to specify the coded names for the fiber and matrix. The format of the data is 
explained in this section so as to enable the user to introduce new contents or to modify existing 
entries as appropriate to his/her needs. Data for four fibers and three matrices are provided in the 
present package. 

The fiber properties are arranged in five physical cards of 80 column length. The first card contains 
a four-character code name of a fiber in format A4. The second to the fifth cards start with a two- 
letter mnemonic to indicate the type of properties that follow. The format on any of these cards is 
A4, 7E10.3, except for the second card. The second card is in the format A3, 16, 7E10.3. The 
mnemonics FP, FE, FT, and FS stand for fiber physical, elastic, thermal, and strength-related 
properties, respectively. The entries on these cards are explained as follows: 
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card 1 four character coded name for fiber 
card 2 FP; Nj, dj, pf 
card 3 FE; Ej^ u E J2 2 , 7/12. 7/23> G f\2> G/23 
card 4 FT; «/22> ^/ii> ^/22> Cf 

card 5 FS; S/y-, Sy C (The remaining entries are open for future modifications.) 

The matrix properties are arranged next after the line OVER END OF FIBER PROPERTIES. The 
properties have essentially the same format as those for fiber property cards. There are, however, six 
physical cards for each matrix material. The mnemonics used are MP, ME, MT, MS, and MV. They 
stand for matrix physical, elastic, thermal, strength-related, and miscellaneous properties, 
respectively. The format for the first card is A4, and the format for the rest of the cards is A3, 
7E10.3. The entries in each card are as follows: 

card 1 four character coded name for matrix 

card 2 MP; p m 

card 3 ME; E m , v m , a m 

card 4 MT ; K m , C m 

card 5 M S ; S in 7-, S mC , t m j, e mC , e mS> t mTOR 
card 6 MV; K„, T gdr 

The data base presently contains properties for T-300 (T300), AS graphite (AS-), S-Glass (SGLA), 
and HMS (HMSF) fibers. The available matrix materials are high-modulus, high-strength (HMHS), 
intermediate-modulus, high-strength (IMHS), and intermediate-modulus, low-strength (IMLS) 
matrices, which are epoxy-type resins. The complete list of properties is shown in appendix C. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 11, 1985 



Appendix A 
List of Code Identifiers 


Engineering Fortran Comment 

symbol symbol code 


^cx 

ACX 

composite axial stiffness; generated in subroutine 
GPCFD2 

A R 

M cx 

RAC 

reduced axial stiffness; computed in subroutine GPCFD2 

BIDE 

Boolean 

true if interply effects are included; input 

C cx 

CPC 

composite coupling stiffness; generated in subroutine 
GPCFD2 

C e] 

RESF 

string with force variables in BLOCK DATA 

C e2 

DISP 

string with displacement variables in BLOCK DATA 

COMSAT 

Boolean 

true if COMSA is executed; input 

CSANB 

Boolean 

true if membrane and bending symmetry exists; input 

Dcx 

FTC 

composite flexural rigidities; generated in subroutine 
GPCFD2 

»§x 

RDC 

reduced bending rigidities; computed in subroutine 
GPCFD2 

D f 

DIAF 

filament equivalent diameter; input 

D v 

DISV, DISVI 

displacement vectors; DISVI is either read in main 
program, or is generated in subroutine COMSA 

Ef> Ecf 

ECF 

filament elastic constants; input 

Efn,n\,m\] 

EFll.ELH.EMll 

filament, ply, and matrix normal moduli; filament and 
matrix moduli input 


EF12,EL12,EM12 

filament, ply, and matrix shear moduli; filament and 
matrix shear moduli input 


ECL 

ply elastic constants; generated in subroutine INHYD 

E cm 

ECM 

matrix elastic constants; generated in subroutine INHYD 

H J 

PL(9,I) 

interply distortion energy coefficient; generated in main 
program 

H kc 

CHK 

array of constituent heat conductivities; input 

h c 

HHC 

composite heat capacity stored in PC(18) and PC(54) 

u 

I,J 

index; generally ply or interply 

IX 

^ cxx,cvy ,cxy 

HK1 1,22,33 

composite two-dimensional heat conductivities in PC(51) 
to PC(53) 

^c\ I,c22,c33 

HK1 1,22,33 

composite three-dimensional heat conductivities along 
the material axes in, PC(15) to PC(17) 

K f.v 

KF, V 

apparent fiber and void volume ratios; input 

^fii.ni.wii 

CHK 

see H kc 

IX 

1X \xx, 1 vy, 1 xy 

XK1,XK2,XK3 

stress concentration factors generated in STRCNF 

k f,m 

KFB,MB 

actual fiber and matrix volume ratios 

k f(, vt 

KFL,VL 

ply apparent fiber and void volume ratios 

k sc 

LSC 

array of limiting conditions; input 

M cx 

MBS 

applied moment; input 

M cT t x 

MSDT 

thermal moments; generated in GPCFD2 

^cM ( x 

MSDH 

hygral moments; generated in GPCFD2 

m 

M 

load condition index 

Hex 

NBS 

applied membrane loads; input 

HcM/X 

NSDH 

hygral force; generated in GPCFD2 

H c T t x 

NSDT 

thermal force; generated in GPCFD2 

Hf 

NFPE 

number of filaments per end; input 


NL 

number of plies; input 
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"k 

NLC 

number of load conditions; input 

N 

1 y ms 

NMS 

number of material systems; input 

"pc 

NPC 

string PROPC length; input 

Npt 

NPL 

string PROP length; input 

NONUDF 

Boolean 

T (true) if Poisson's ratio difference chart is to be 
suppressed 

Pc 

PC 

composite properties array; generated in GACD3 and 
GPCFD2 

Pep 

PROPC 

string PROPC; composite property identifiers in 
GDCFD2 

p f 

PL 

ply property array; portions generated in ail parts of the 
program 

Pfp 

PROP 

string PROP; ply properties identifiers in main program 

Qf,i,p,r,s 

QF,I,P,R,S 

indices to print out string PROP 

R 

R 

transformation matrix; GACD3, GPCFD2, and COMSA 

RINDV 

Boolean 

T (true) if displacements are read in; input 

Sg nr* et c. 

PL(51) to PL(59J) 

ply limit stresses; generated in GLLSC 

U 

TL 

ply thickness; input 

w cb 

w xx 

composite local curvatures relative to the structural axes 

a c 

CTE 

composite coefficient of thermal expansion; three- 
dimensional in PC(12) to PC(14), two-dimensional in 
PC(48) to PC(50) 

a f, a, m 

VAF,AL,AM 

filament, ply, and matrix thermal coefficients of 
expansion; input 

0e y i>e 

VCF 

correlation factors for ply thermoelastic properties and 
strain magnification factors; set to unity in COMSA 

h 

BTA 

correlation factors for ply heat conductivity; set to unity 
in COMSA 

Ps 

BET 

correlation factors for ply strength; set to unity in 
COMSA 

6 ( 

PL(8,I) 

interply layer thickness; generated in INHYD 

^csx 

UX 

reference plane membrane strain; solved in terms of N cx 
or input 


EPS,PL(74) to PL(66,I) 

ply strains; generated in COMSA 

Ocs 

THCS 

angle between composite material and structural axes; 
input 


THLC 

angle between ply material and composite axes; input 

v f\2 y fl2, m\2 

NUF12,L12,M12 

filament, ply, and matrix Poisson’s ratios; input 

7T 

PIE 

constant; input 

Pf,mJ 

RHOF,M,L 

filament and matrix weight density; input and generated 
in FIBMT, HTM, and COMPP 

°f> °rn 

SF, SM 

microstresses in fibers and matrices generated in 
MCRSTR 

ti 

XPL,XSL,YPL,YSL 

boundary zone decay length; generated in the main 
program and paired to EDGSTR 


SIGL,PL(67) to PL(69,I) 

ply stress; generated in COMSA 



Appendix B 
Sample Input/Output 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Item 1 

I C A N 


******** 
******** 
** 

* * 

** 

** 

** 


******************* ****************************************** 
************************************************************* 


** 

* * 

* * 

* * 

** 

******** 

******** 


NASA LERC. CLEVELAND 
SEPTEMBER 1983 


**** 
**** 
** 
** 
* * 
** 
** 


X* 

SSSSSSSSSS 

CCCCCCCCC 

AAAAAAAA 

NNN 


NNN 

XX 

XX 

SSSSS5S5SSS5 

CCCCCCCCCCC 

AAA AAA AAA A 

NNNN 


NNN 

XX 

** 

SSSS 

CCC 

AAA 

AAA 

NNMNN 

NNN 

XX 

** 

SSSS 

CCC 

AAA 

AAA 

NNNNNN 

NNN 

X* 

* * 

S5SS 

CCC 

AAA 

AAA 

NNN NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAAAAAAAAAAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA AAA AAA AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

X* 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNN 

NNN 

XX 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNNNNN 

X* 

XX 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNNNNN 

X* 

X* 

SSSS 

CCC 

AAA 

AAA 

NNN 

NNMNN 

X* 

X* 

S5SSS5SSSSSS 

CCCCCCCCCCC 

AAA 

AAA 

NNN 


NNNN 

XX 

XX 

XX 

SSSSSSSSSS 

CCCCCCCCC 

AAA 

AAA 

NNN 


NNN 

XX 

XX 

XX 

XX 

INTEGRATED 

COMPOSITES ANALYZER 

BY P . 

, L .N. 

MURTHY 

XX 

XX 


************************************************************* 

************************************************************* 


** 
* * 
** 
xx 
xx 
XXX* 
xxxx 


Item 2 

I C A N: COORDINATE SYSTEMS 


/ /I 

/ /A | 

/ /AB t 

/ /ABB | 

/. /ABBA/ 

/. . /ABBA/ 

/. . . THETA /ABBA/ 

/. . . . /ABBA/ 

/ /ABBA/ 

ABBA/ 

1000000000000000000000000000000 i BB A/ 

| XXXXXXXXXXXXXXXXXXXXXXXXXXXXXX i BA/ 

| XXXXXXXXXXXXXXXXXXXXXXXXXXXXXX | A/ 

i 0 00000000 000000000000000000000 I / 


/ 


/I 

/B/- 
/B/ 
/B/ 

Y / /B/ 

/ /B/ 

/ /B/ 

/ /B/ 

/ /B/ 

B/ 

i ooooooooooooooooooooooooooooool/ 


LAMINATE STRUCTURAL AXES 


PLY MATERIAL AXES 



Item 3 


I C A H 

INPUT DATA ECHO 


FOUR PLY SYMMETRIC LAMINATE. ICAN SAMPLE INPUT DATA. 


STDATA 
T 
F 
F 
F 
T 


COMSAT 

CSANB 

BIDE 

RINDV 

NONUDF 


PLY 1 


1 

70.00 

70.0 

. 0 

0.0 

.010 

PLY 2 


2 

70.00 

70.0 

. 0 

90 . 0 

.005 

PLY 3 


2 

70.00 

70.0 

. 0 

90.0 

.005 

PLY A 


1 

70.00 

70.0 

. 0 

0.0 

.010 

MATCRDAS — IMLS 

.55 


. 02 

AS— IMLS 

0.0 

.57 

.03 

MATCRDSGL AHMHS 

.55 


.01 

AS— IMHS 

0.9 

.57 

.01 

PLOAD 1000. 

0 . 0 


0.0 

0.0 


NX , NY > 1 

NXY, THCS 

PLOAD 0.0 

0 . 0 


0 . 0 



MX P MY ,1 

MXY 

PLOAD 0.0 

0 . 0 





DMX/QX 

, DMY/QY , PR5S 


Item 4 

SUMMARY OF INPUT DATA 

FOUR PLY SYMMETRIC LAMINATE. ICAN SAMPLE INPUT DATA. 


- - - CASE CONTROL DECK - ■ 
NUMBER OF LAYERS NL = A 
NUMBER OF LOADING CONDITIONS NLC = 1 
NUMBER OF MATERIAL SYSTEMS NMS = 2 


COMSAT CSANB BIDE 
T F F 


RINDV NONUDF 
F T 


LAMINATE CONFIGURATION 


PLY 

NO 

MID 

DELTA T 

DELTAM 

THETA 

T-NESS 

PLY 

1 

1 

0.000 

0.0% 

0.0 

0 . 010 

PLY 

2 

2 

0.000 

0.0% 

90.0 

0 . 005 

PLY 

3 

2 

0.000 

0.0% 

90.0 

0.005 

PLY 

A 

1 

0 . 000 

0 . 0% 

0.0 

0.010 


- - - 

COMPOSITE MATERIAL SYSTEMS 

. - - 




MATCRD 

MID 

PRIMARY 

VFP 

VVP 

SECONDARY 

V5C 

VFS 

vvs 

MATCRD 

1 

AS— IMLS 

0.55 

0.02 

AS — IMLS 

0.00 

0.57 

0.03 

MATCRD 

2 

SGLAHMHS 

0.55 

0.01 

AS — IMHS 

0 .90 

0.57 

0.01 


LOADING CONDITIONS - - - 


PRESCRIBED LOADS FOR 

THE LOAD 

CONDITION 1 


INPLANE LOADS 

NX 

= 

1000.0000 

LB/IN 


NY 

z 

0 . 0000 

LB/IN 


NXY 

z 

0 . 0000 

LB/IN 

3ENDING LOADS 

MX 

= 

0 . 0000 

LB. IN/IN 


MY 

= 

0 . 0 000 

LB . IN/IN 


MXY 

= 

0 . 0000 

LB . IN/IN 

TRANSVERSE LOADS 

DMX/QX 

z 

0.0000 

LB/IN 


DMY/QY 

- 

O .0000 

LB/IN 

TRANSVERSE PRESSURE 

PU 

- 

0.0000 

LB/SQ. IN 

TRANSVERSE PRESSURE 

PL 

z 

0 .0000 

LB/SQ. IN 
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Item 5(a) 


--> CONSTITUENT PROPERTIES: ECHO FROM DATA BANK. <— 
PRIMARY FIBER PROPERTIES; AS— FIBER 


ORIGINAL 

OF POOR 


PAGE is 

Quality 



PRIMARY MATRIX PROPERTIES; IMLS MATRIX. DP.Y RT . PROPERTIES. 


1 

ELASTIC MODULUS 

EMP 

0 .5000E 06 

2 

SHEAR MODULUS 

GMP 

0 . 1773E 06 

3 

POISSON* S RATIO 

HUMP 

0 . 4100E 00 

4 

THERM. EXP. COEF. 

CTSMP 

0 . 57 0 0E- 04 

5 

DENSITY 

RHOMP 

0 . 4600E-01 

6 

HEAT CAPACITY 

CMPC 

0.2500E 00 

7 

HEAT CONDUCTIVITY 

KMP 

0.1250E 01 

8 

STRENGTHS 

SMPT 

0 .7000E 04 

9 


SMPC 

0 . 2100E 05 

10 


SMPS 

0 . 7000E 04 

11 

MOISTURE COEF 

BTAMP 

0 . 40 00E-02 

12 

DIFFUSIVITY 

DIFMP 

0 . 2000E-03 


Item 5(b) 

PRIMARY COMPOSITE PROPERTIES; 55/ 43 AS — /IMLS 

BASED ON MICROMECHANICS OF INTRAPLY HYBRID COMPOSITES: ELASTIC AND THERMAL PROPERTIES. 


FIBER VOLUME RATIO - 0.550 MATRIX VOLUME RATIO - 0.430 VOID VOLUME RATIO - 0.020 

VOID CONDUCTIVITY - 0.22499990E 00 


1 

ELASTIC MODULI 

EPC1 

0.1726E 08 

2 


EPC2 

Q.1127E 07 

3 


EPC3 

0.1127E 07 

4 

SHEAR MODULI 

GPC12 

0.5470E 06 

5 


GPC23 

0.3238E 06 

6 


GPC13 

0 . 547 OE 06 

7 

POISSON’S RATIO 

NUPC12 

0.2945E 00 

8 


NUPC23 

0 . 482 IE 00 

9 


NUPC13 

0 . 2945E 00 

10 

THERM. EXP. COEF. 

CTEPC1 

0 . 1418E-06 

11 


CTEPC2 

0 .2464E-04 

12 


CTEPC3 

0 .2464E-04 

13 

DENSITY 

RHOPC 

0 .5443E-01 

14 

HEAT CAPACITY 

CPC 

0 . 1991E 00 

15 

HEAT CONDUCTIVITY 

KPC1 

0 . 3195E 03 

16 


KPC2 

0 . 3702E 01 

17 


KPC3 

0 . 3702E 01 

18 

STRENGTHS 

SPC1T 

Q.2228E 06 

1 9 


SFC1C 

0.8764E 05 

20 


SPC2T 

0 . 5006E 04 

21 


SPC2C 

0.1502E 05 

22 


SPC12 

0.5126E 04 

23 

MOIST. DIFFUSIVITY 

DPC1 

Q .8600E-04 

24 


DPC2 

0 .5168E-04 

25 


DPC3 

0 .516SE-04 

26 

MOIST. EXP. COEF. 

BTAPC1 

0 . 4981E-04 

27 


BTAPC2 

0 . 1452E-02 

23 


BTAFC3 

0 . 1452E-02 

29 

FLEXURAL MODULI 

EPC1F 

0.1726E 08 

30 


EPC2F 

0 . 1127E 07 

31 

STRENGTHS 

SFC23 

0.3983E 04 

32 


SPC1F 

0.1572E 06 

33 


SPC2F 

0.9387E 04 

34 


SPCSB 

0.7689E 04 

35 

PLY THICKNESS 

TPC 

0 . 5 0 00E-02 

36 

INTERPLY THICKNESS 

PLPC 

0 .585QE-04 

37 

INTERFIBER SPACING 

PLPCS 

0 .5850E-04 
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Item 5(c) 

— > CONSTITUENT PROPERTIES: ECHO FROM DATA BANK. < — 
PRIMARY FIBER PROPERTIES; SOLA FIBER 


1 

ELASTIC MODULI 

EFP 1 

0.1240E 08 

2 


EFP2 

0.1240E 08 

3 

SHEAR MODULI 

GFP12 

0.5170E 07 

4 


GFP23 

0.5170E 07 

5 

POISSON * 5 RATIO 

NUFP12 

0.2000E 00 

6 


NUFP23 

0.2000E 00 

7 

THERM. EXP. COEF. 

CT EFP 1 

0.2800E-05 

8 


CTEFP2 

0.28COE-05 

9 

DENSITY 

RHOFP 

0 . 90 (JOE- 01 

10 

NO. OF FIBERS/END 

NFP 

0.2040E 03 

11 

FIBER DIAMETER 

DIFP 

0 . 3600E-03 

12 

HEAT CAPACITY 

CFPC 

0 . 1700 E 00 

13 

HEAT CONDUCTIVITY 

KFP1 

0.7500E 01 

14 


KFP2 

0.7500E 01 

15 


KFP 3 

0.7500E 01 

16 

STRENGTHS 

SFPT 

0.3600E 06 

17 


SFPC 

0.3000E 06 

PRIMARY 

MATRIX PROPERTIES; 

HMHS MATRIX. 

DRY RT. PROPERTIES 

1 

ELASTIC MODULUS 

EMP 

0.7500E 06 

2 

SHEAR MODULUS 

GMP 

0.2778E 06 

3 

POISSON’S RATIO 

NUMP 

0 . 350 OE 00 

4 

THERM. EXP. COEF. 

CTEMP 

0.400 0 E-04 

5 

DENSITY 

RHOMP 

0 .4500E-01 

6 

HEAT CAPACITY 

CMPC 

0 . 2500 E 00 

7 

HEAT CONDUCTIVITY 

KMP 

0.1250E 01 

8 

STRENGTHS 

SMPT 

0.2000E 05 

9 


SMPC 

0.5000E 05 

1 0 


SMPS 

0.1500E 05 

11 

MOISTURE COEF 

BT AMP 

0 . 4000E-02 

12 

DIFFUSIVITY 

DIFMP 

0.2000E-03 



Item 5(d) 


PRIMARY 

COMPOSITE PROPERTIES; 

55/ 44 SGLA/HMHS 


BASED ON MICROMECHANICS OF INTRAPL Y HYBRID COMPOSITES: ELASTIC AND THERMAL PROPERTIES. 

FIBER VOLUME RATIO - 0.550 MATRIX VOLUME RATIO - 0.440 VOID VOLUME RATIO - 0 010 
VOID CONDUCTIVITY - 0.22499990E 00 


1 

ELASTIC MODULI 

EPC1 

0.7150E 07 

2 


EPC2 

0.2473E 07 

3 


EPC3 

0.2473E 07 

4 

SHEAR MODULI 

GPC12 

0.9314E 06 

5 


GPC23 

0.5792E 06 

6 


GPC13 

0.9314E 06 

7 

POISSON’S RATIO 

NUPC12 

0.2675E 00 

8 


NUPC23 

0.3778E 00 

9 


NUPC13 

0.2675E 00 

10 

THERM. EXP. COEF. 

CTEPC1 

0 . 4488E-05 

11 


CTEPC2 

0.1580E-04 

12 


CTEPC3 

0.1580E-04 

13 

DENSITY 

RHOPC 

0. 6930 E- 01 

14 

HEAT CAPACITY 

CPC 

0.1929E 00 

15 

HEAT CONDUCTIVITY 

KPC1 

0.4675E 01 

16 


KPC2 

0.2750E 01 

17 


KPC3 

0.2750E 01 

18 

STRENGTHS 

SPC1T 

0.2076E 06 

19 


SPC1C 

0.1730E 06 

20 


SPC2T 

0.1256E 05 

21 


SPC2C 

0.3140E 05 

22 


5PC12 

0.1047E 05 

23 

MOIST. DIFFUSIVITY 

DPC1 

0.8800E-04 

24 


DPC2 

0 .5168E-Q4 

25 


DPC3 

0.5168E-04 

26 

MOIST. EXP. COEF. 

BTAPC1 

0 . 1846E-03 

27 


BTAPC2 

0. 1379E-02 

28 


BTAPC3 

0. 1379E-02 

29 

FLEXURAL MODULI 

EPC1F 

0.7150E 07 

30 


EPC2F 

0.2473E 07 

31 

STRENGTHS 

SPC23 

0 . 651 OE 04 

32 


SPC1F 

0.2359E 06 

33 


SPC2F 

0.2243E 05 

34 


SPCSB 

0 . 157 OE 05 

35 

PLY THICKNESS 

TPC 

0. 5000 E- 02 

36 

INTERPLY THICKNESS 

PL PC 

0 .7020E-04 

37 

INTERFIBER SPACING 

PL PCS 

0. 7020 E- 04 
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Item 5(e) 


--> CONSTITUENT PROPERTIES: ECHO FROM DATA BANK. ■>- 
SECONDARY FIBER PROPERTIES; AS— FIBER 


1 

ELASTIC MODULI 

EFS1 

0 . 310CE 03 

2 


EFS2 

0.2000E 07 

3 

SHEAR MODULI 

GFS12 

0.2000E 07 

4 


GFS23 

0 . 1000E 07 

5 

POISSON’S RATIO 

NUFS12 

0.2000E 00 

6 


NUFS23 

0.2500E 00 

7 

THERM. EXP. COEF. 

CTEFS1 

-0.5500E-06 

8 


CTEFS2 

0 . 5600E-05 

9 

DENSITY 

RHOFS 

0 . 6300E-01 

10 

NO. OF FIBERS/END 

NFS 

0.1000E 05 

11 

FIBER DIAMETER 

DIFS 

0 . 3000E-03 

12 

HEAT CAPACITY 

CFSC 

0 . 1700E OO 

13 

HEAT CONDUCTIVITY 

KFS1 

0.5800E 03 

14 


KFS2 

0 . 5800E 02 

15 


KFS3 

0 . 58 0 0 E 02 

16 

STRENGTHS 

SFST 

0 . 4000E 06 

17 


SFSC 

0.4000E 06 


SECONDARY 

MATRIX PROPERTIES; 

IMHS MATRIX 

. DRY RT. } 

1 

ELASTIC MODULUS 

EMS 

0 . 5000E 06 

2 

SHEAR MODULUS 

GMS 

0.1852E 06 

3 

FOISSON’S RATIO 

HUMS 

0 . 35 0 0 E 00 

4 

THERM. EXP. COEF. 

CTEMS 

0 . 3600E-04 

5 

DENSITY 

RHOMS 

0 .4400E-01 

6 

HEAT CAPACITY 

CtlSC 

0.2500E 00 

7 

HEAT CONDUCTIVITY 

KMS 

0.1250E 01 

8 

STRENGTHS 

snsT 

0 . 1500E 05 

9 


SMSC 

0.3500E 05 

1 0 


SMSS 

0 . 1300E 05 

11 

MOISTURE COEF 

BTAtlS 

0 . 40 0 0E-02 

12 

DIFFUSIVIT Y 

DIFMS 

0 . 200 0E-03 



Item 5(f) 

SECONDARY 

COMPOSITE PROPERTIES; 

57/ 42 AS- 

-/IMHS 


BASED ON MICROMECHANICS OF INTRAPLY HYBRID COMPOSITES: ELASTIC AND THERMAL 


FIBER VOLUME RATIO - 0.570 MATRIX VOLUME RATIO - 0.420 VOID VOLUME 

VOID CONDUCTIVITY - 0.22499990E 00 


1 

ELASTIC MODULI 

ESC1 

0.1788E 08 

2 


ESC2 

0.1153E 07 

3 


ESC3 

0 . 1153E 07 

4 

SHEAR MODULI 

GSC12 

0 . 588 0E 06 

5 


GSC23 

0.3458E 06 

6 


GSC13 

0 . 588 0E 06 

7 

POISSON'S RATIO 

NUSC12 

0.2645E 00 

8 


NUSC23 

0 . 4294E 00 

9 


NUSC13 

0 . 2645E 00 

10 

THERM. EXP. COEF. 

CTESC1 

-0 . 1280E-06 

11 


CTESC2 

0 . 1 605E-04 

12 


CTESC3 

0 . 1605E-04 

13 

DENSITY 

RHOSC 

0 . 5439E-01 

14 

HEAT CAPACITY 

CSC 

0 . 1 972E 00 

15 

HEAT CONDUCTIVITY 

KSC1 

0.3311E 03 

16 


KSC2 

Q.3918E 01 

17 


KSC3 

0.3918E 01 

18 

STRENGTHS 

SSC1T 

0.2307E 06 

19 


SSC1C 

0.1568E 06 

20 


SSC2T 

0.1026E 05 

21 


SSC2C 

0 .2394E 05 

22 


SSC12 

0.9369E 04 

23 

MOIST. DIFFUSIVITY 

DSC1 

0 .8400E-04 

24 


DSC2 

0 .4900E-04 

25 


DSC3 

0 .4910E-04 

26 

MOIST. EXP. COEF. 

BTASC1 

0 .46V8E-04 

27 


BTASC2 

0 . 1319E-02 

28 


3TASC3 

0 . 1319E-02 

29 

FLEXURAL MODULI 

ESC IF 

0.1788E 03 

30 


ESC2F 

0.1153E 07 

31 

STRENGTHS 

SSC23 

0.7424E 04 

32 


SSC1F 

0.2334E 06 

33 


SSC2F 

0.1796E 05 

34 


SSCSB 

0 . 1405E 05 

35 

PLY THICKNESS 

TSC 

0 .5000E-02 

36 

INTERPLY THICKNESS 

PLSC 

0 . 5215E-04 

37 

INTERFIBER SPACING 

PLSCS 

0 . 5215E-04 


PROPERTIES. 
RATIO - 0.010 


59 



Item 5(g) 

HYBRID COMPOSITE PROPERTIES; 60/40 SGLA/HMHS/flS — /IMHS 

BASED ON MICROMECHANICS OF INTRAPLY HYBRID COMPOSITES: ELASTIC AND THERMAL PROPERTIES. 


PRIMARY COMPOSITE VOLUME RATIO - 0.600 SECONDARY COMPOSITE VOLUME RATIO - 0.400 


1 

ELASTIC MODULI 

EHC1 

0.1144E 08 

2 


EHC2 

0 . 1696E 07 

3 


EHC3 

0.1945E 07 

4 

SHEAR MODULI 

GHC12 

0.7551E 06 

5 


GHC23 

0 . 456 IE 06 

6 


GHC13 

0.7941E 06 

7 

POISSON'S RATIO 

NUHC12 

0.2663E 00 

8 


NUHC23 

0.3985E 00 

9 


NUHC13 

0.2689E 00 

10 

THERM. EXP. COEF. 

CTEHC1 

0.1603E-05 

11 


CTEHC2 

0 . 1601E-04 

12 


CTEHC3 

0 . 1634E-04 

13 

DENSITY 

RHOHC 

0 . 6334E-01 

14 

HEAT CAPACITY 

CHC 

0.1943E 00 

15 

HEAT CONDUCTIVITY 

KHC1 

0.1352E 03 

1 6 


KHC2 

0.2305E 01 

17 


KHC3 

0.2305E 01 

18 

STRENGTHS 

SHC1T 

0.2168E 06 

1 9 


SHC1C 

0.1665E 06 

20 


SHC2T 

0.9915E 04 

21 


SHC2C 

0.2314E 05 

22 


SHC12 

0.1195E 05 

23 

MOIST. DIFFUSIVITY 

DHC1 

0 .8736E-04 

24 


DHC2 

Q.5117E-04 

25 


DPC3 

0 .5117E-04 

26 

MOIST. EXP. COEF. 

BTAHC1 

0 . 9858E-04 

27 


BTAHC2 

0 .8565E-03 

28 


BTAHC3 

0 . 1455E-02 

29 

FLEXURAL MODULI 

EHC1F 

0.1144E 08 

30 


EHC2F 

0.1696E 07 

3 1 

STRENGTHS 

SHC23 

0.1019E 05 

32 


SHC1F 

0.2355E 06 

33 


SHC2F 

0.1735E 05 

34 


SHCSB 

0.1793E 05 

35 

PLY THICKNESS 

THC 

0 . 5000E-02 

36 

INTERPLY THICKNESS 

PLHC 

0.5215E-04 

37 

INTEP.FIBER SPACING 

PLHCS 

0.5215E-04 

3 8 

FIBER VOL. RATIO 

VFH 

0.5580E 00 

39 

MOISTURE CONTENT 

M 

0.0000 

40 

MATRIX VOL. RATIO 

VMH 

0.4320E 00 


60 



ORIGINAL PAGE IS 
OF POOR QUALITY 


Item 6 


3-D COMPOSITE STRAIN STRESS TEMPERATURE MOISTURE RELATIONS - STRUCTURAL AXES 


-1- 

-z- 

-3- 

-4- 

-5- 

-6- 

-DT- 


0 . 6976E-07 

-0 . 5 9 52 E- 0 8 

-0 .2727E-07 

0.0 0 0 0 

0 . 0 0 0 0 

0 . 3255E-1 3 

0 . 1102E-05 

0 . 

-0.5952E-08 

0 . 19626-06 

-0.8485E-07 

0.0000 

0.0000 

-0.I464E-11 

0 . 5805E-05 

0 

-0 . 2727E-07 

-0 .8485E-07 

0.5614E-06 

0.0000 

0.0000 

0 .6682E-12 

0 .2839E-04 

0 

0.0000 

0.0 0 0 0 

0 .0000 

0 .2139E-05 

0.5229E-12 

0 .0000 

0 .0 0 00 


0. 0000 

0.0000 

0.0 0 0 0 

0.5229E-12 

0 . 1 935E-05 

0 . 0000 

0.0000 


0 .3255E-I3 

-0 . 1464E-11 

0 . 6682E-12 

0 .0000 

0 .0000 

0. 1622E-05 

-0 .4330E-10 

-0 


3-D COMPOSITE STRESS STRAIN RELATIONS - STRUCTURAL AXES 




-1- 

-2- 

-3- 

-4- 

-5- 

-6- 


1 

0.1473E 08 

0 . 8093E 06 

0 . 8377 E 06 

0.0 0 0 0 

0.0000 

0 .8960E-01 


2 

0.8093E 06 

0.5499E 07 

0 . 87 0 3E 06 

0.0000 

0.000 0 

0.4587E 01 


3 

0.8377E 06 

0 .S703E 06 

0 . 1953E 07 

0.0000 

0.0000 

0 .3603E-01 


4 

0.0000 

0.0000 

0 .0000 

0.4676E 06 

0.I263E 00 

0.0000 


5 

0.0000 

0 . 0000 

0.0000 

0.1263E 00 

0.5167E 06 

0.0000 


6 

0 .896 0E-01 

0 .4587E 01 

■0 . 3603E-01 

0.0000 

0 . 0000 

0.6I64E 06 



MAT 9 CARD FOR MSC/NASTR AN SOLID ELEMENTS 



G56 i G66 

0 . 0 0 00 00 00 
0 . 00000000 


0.54987690E 07 
0.00000000 


0 . 87 0 32 4 0 6 E 
0 . 61636813E 


Item 7 

COMPOSITE PROPERTIES 


COMPOSITE P ROPERTI 
LINES 1 TO 31 3-D 
LINES 33 TO 62 2 ~ 


ES - VALID ONLY FOR CONSTANT TEMPERATURE THROUGH THICKNESS 
COMPOSITE PROPERTIES ABOUT MATERIAL AXES 
D COMPOSITE PROPERTIES ABOUT 5TRUCTURAl_AXES 


1 

RHOC 

0 . 5740E-01 

2 

TC 

0 . 3000E-01 

3 

ecu 

0.1473E 08 

4 

CC12 

0 .8093E 06 

5 

CC13 

0 .8377E 06 

6 

CC22 

0 . 54 9 9 E 07 

7 

CC2 3 

0 . 87 0 3 E 06 

8 

CC33 

0 . 1953E 07 

9 

CC44 

0 . A 6 7 6 E 06 

10 

CCS 5 

0 .5167E 06 

u 

CC6 6 

0 . 6164E 06 

12 

CTE11 

0 . 1102E-05 

13 

CTE22 

0 . 5805E-05 

14 

CTE33 

0 .2839E-04 

15 

HK1 1 

0.2138E 03 

16 

HK22 

0.4755E 02 

1 7 

HK33 

0 .3083E 01 

18 

HHC 

0 . 1 97 5E 00 

19 

EC1 1 

0 . 1433 E 08 

20 

EC22 

0 . 5 0 98 E 07 

21 

EC 3 3 

0 . 1 7 3 1 E 07 

22 

EC23 

0 . 467 6 E 06 

23 

EC3 1 

0 . 5167E 06 

26 

EC12 

0 . 6I64E 06 

25 

HUC12 

0 . S531E-01 

26 

NUC21 

0 . 30 36 E- 0 1 

27 

H'JC 1 3 

0 . 390SE 00 

28 

NUC31 

0 .4856E-01 

29 

NUC23 

0 . 4 326 E 00 

30 

NUC32 

0 . 1511E 00 

31 

ZCGC 

0 . 1500E-01 


32 

32DEC 

0 .0000 

33 

CC11 

0.1215E 08 

34 

CC12 

0 . 3746 E 06 

35 

CC13 

0.1067E 00 

36 

CC22 

0.4610E 07 

37 

CC23 

0.4055E 01 

38 

CC33 

0.6164E 06 

39 

ECU 

Q.1212E 08 

60 

EC22 

0.4598E 07 

41 

EC 1 2 

0.6164E 06 

42 

NUC12 

0. 8126 E- 01 

43 

NUC21 

0. 3084E-01 

44 

CSN13 

-0 . 3615E-G6 

65 

CSN31 

-0 . 1839 E- 0 7 

66 

CSN23 

0 . 6574E-05 

67 

CSN32 

0 . 8811E-Q6 

48 

CTE11 

0 . 1102E-05 

49 

CTE22 

0. 5805 E- 05 

50 

CTE12 

-0. 4330 E- 10 

51 

HK11 

0 .2138E 03 

52 

HK22 

0.4755E 02 

53 

HK12 

-0 . 5785E-04 

54 

HHC 

0 . 1975E 00 

55 

DPC11 

0 . 3573E-02 

56 

DPC22 

0 , 3578E-02 

57 

DPC33 

0.2000E 05 

58 

DPC12 

3.0000 

59 

BTAC11 

0 . 100 9E-03 

60 

BTAC22 

0 . 337 0 E-03 

61 

BTAC33 

0 . 18 59 E- 02 

62 

BTAC12 

-0 . 2425E-08 


-DM- 

1009E-03 
3370E-03 
. I859E-02 
0 . 0000 
0 .0000 
.2425E-0S 


06 

06 


61 



Item 8 


FORCES 



FORCE DISPLACEMENT RELATIONS 



DI5PL 

T-F0RCE5 

H-F0RCE5 

NX 

0 . 3644 E . 06 

Q.1124E 05 

0.3201 E— 0 2 

0 . 0000 

0.0000 

-0 .7276E-11 

UX 

0.0000 

0 .0 0 0 0 

NY 

0.1124E 05 

0.1383E 06 

0.1217E 00 

0 . 0000 

-0.1373E-03 

-0.2619E-09 

VY 

0.0000 

0 . 0 0 0 0 

NXY 

0 . 3201E-02 

Q.1217E 00 

0 . 184 9 E 05 

-0 . 7276E-11 

-0 .2619E-09 

0.0000 

VXPUY 

0.0000 

0 . 0 0 0 0 

MX 

0 .0000 

0 . CCCO 

- 0 . 7 27 6 E- 1 1 

0.3776E 02 

0.7610E 00 

0. 2667 E- 07 

wxx 

0.0000 

0 .0000 

MY 

0 .0000 

-0.1373E-03 

-0 .2619E-09 

0.7610E 00 

0.3419E 01 

0. 1014E-05 

WYY 

0 .0000 

0.0000 

MXY 

-0 .7276E-11 

-0 . 2619E-09 

0 . 0 0 0 0 

0 .2667E-07 

0 .1014E-05 

0.1248E 01 

WXY 

0.0000 

0 .0000 


Item 9 

REDUCED STIFFNESS MATRIX REDUCED BENDING REGIDITIES 


Q.36441E 06 
0 . 11238c 05 
0 . 32306E-02 

0. 1 1238E 
0. 13S30E 
0. 12165E 

05 0.32006E-02 

06 0 . 12165E 00 

00 0 . 184 9 1 E 05 

0.37763E 02 
0.76104E 00 
0 .26672E-07 

0.76104E 00 
0 . 34186 E 01 
0 . 1Q13SE-05 

0.26672E-07 
0 . 10138E-05 
Q.12481E 01 




Item 10 



C M E US 

E F U L 

DATA FOR F.E. 

ANALYSIS 




COMPOSITE THICKNESS FOR F.E. ANALYSIS = 0.30000E-Q1 

PROPERTIES FOR F.E. ANALYSIS E 1 1 , El 2 . El 3 , E22 , E23 . E33 PROPERTIES SCALED BY 10**-6 
3 . 32533 E- 01 - 0 . 6 7 0 6 9 E- 0 2 -0.29839E-Q7 0.21747E 0 0 0.14296E-05 0.16224E 01 


BENDING EQUIVALENT PROPERTIES NUCXY, NCYX, ECXX, ECYY , GCXY 
3.22261E 00 0.20153E-0I 0.1670SE 08 0.15126E 07 0.55473E 06 

N A £ TRAN MEMBRANE EQUIVALENT ELASTIC COEFFICIENTS 011,012,013,022,023,033 
3.12147E 08 0 . 37 <46 2E 06 0.10669E 00 0.46Q99E 07 0.40551E 01 0.61637E 06 

NASTRAN BENDING EQUIVALENT ELASTIC COEFFICIENTS 01 1 , G12 , 013 , G22 . G23 , G33 
3.167G4E 08 0.33824E 06 0.11854E-01 0.15194E 07 Q.45056E 00 0.55473E 06 


Item 11 

DI5P. DISPLACEMENT FORCE RELATIONS COMBINED FORCES 


1 

0 . 27 51 E- 02 

-1- 

0.2751E-G5 

-2- 

-0 .2236E-Q6 

-3- 

0.9946E-12 

-4- 

0.1818E-12 

-5- 

-O.9021E-11 

-6- 

-0 . 2356E-16 

0.1000E 04 

2 

-0 .2236E-03 

-0 .2236E-06 

0 .7249E-05 

-0 .4765E-10 

-0 . 5895E-1 1 

0.2925E-09 

0 . 1283E-14 

0.0000 

3 

0. 9946E-09 

0 . 9946E-12 

-0 .4765E-10 

0 -5408E-04 

-0 . 3466E-16 

0 . 2237E-14 

-0.1181E-19 

0.000 0 

4 

0 . 1818E-09 

0 . 1818E-12 

-0 .5895E-11 

-0 .3466E-16 

0.2660E-01 

-0.5922E-02 

0.4241E-08 

0.0000 

5 

-0 . 9 0 21 E- 08 

-0 . 9021E-11 

0.2925E-09 

0.2237E-14 

-0 .5922E-02 

0.2938E 00 

-0 .2385E-06 

0 .0000 

6 

-0 .2356E-13 

-0. 2356E-16 

0. 1283E-14 

-0. 1181E-19 

0 . 4241 E-08 

-0.2385E-06 

0.8012E 00 

0.0000 


NOTE: THE DISPLACEMENTS ARE REFERENCE PLANE MEMBRANE STRAINS (UX , VY , VXPUY) AND CURVATURES CUXX , WYY , WXY) 
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Item 12 


PLY hygrqthermomechanical properties/response 


FOR LOAD CONDITIONS 
MEMBRANE LOADS NBS( X , Y , XY~M) ARE 

BENDING LOADS MBS( X , Y , XY-M) ARE 

QXZrQYZ AND APPLIED PRESSURES ARE 
NOTE : NO MOISTURE OR TEMPERATURE 


1 0 0 0 . 
0 . 
0 . 


0 . 0 . 

0 . o . 

0 . 0 . 0. 


LAYER PROPERTIES , RONS-PROPERTY , COLUMNS-LAYER 


PLY NUMBER 1 

MATERIAL SYSTEM AS--/IMLS 
AS — /IMLS 

ORIENTATION 0 . 0 


2 

SGLA/HMHS 
AS--/IMHS 
90 . 0 


3 

SGLA/HMHS 
AS--/IMHS 
90 . 0 


A 

AS — /IMLS 
AS- -/IMLS 
0 . 0 


1 

KV 

2 

KF 

3 

KFB 

4 

KM 

5 

KMB 

6 

P.HOL 

7 

TL 

3 

DELTA 

9 

ILDC 

10 

ZB 

11 

ZGC 

12 

THCS 

13 

THLC 

14 

THLS 

15 

SCI 1 

16 

SC12 

17 

SCI 3 

18 

SC22 

19 

SC23 

20 

SC3 3 

21 

SC44 

22 

SC55 

23 

SC6 6 

24 

CTE1 1 

25 

CTE22 

26 

CTE33 

27 

HK11 

28 

HK22 

"29 

HK33 

30 

HCL 

31 

ELI 1 

32 

EL22 

33 

EL 3 3 

34 

GL23 

35 

GL1 3 

36 

GL12 

137 

NUL12 

38 

NUL21 

39 

NUL13 

40 

NUL31 


0.2000E-01 
0 . 55 0 0 E 00 
0.5390E 00 
0 . 4500E 00 
0 . 4410E 00 
0 . 5443E-01 
0 . 10Q0E-O1 
0 . 5850E-Q4 
0.0000 
0 .5000E-02 
-0 . 1000E-01 
0 .0 0 0 0 
0.0000 
0.0000 
0.2101E 08 
0.7797E 06 
0.7797E 06 
0.1631E 07 
0.8713E 06 
0 . 1776E 07 
0.3238E 06 
0 . 547 0E 06 
0 . 54 7 0 E 06 
0 . 1413E-06 
0 . 2464E-04 
0 . 2464E-Q4 
0 . 3195E 03 
0 . 37 0 2E 01 
0 . 3702E 01 
0 . 1991E 00 
0 . 1 7 26E 08 
0 . 1127E 07 
0 . 1127E 07 
0 . 3238E 06 
0 . 547 0E 06 
0 . 54 7 0 E 06 
0.2945E 00 
0 . 1922E-01 
0.2945E 00 
0 . 1922E-01 


0 . 1000E-01 
0 . 558 0E 00 
0.5524E 00 
0.4420E 00 
0 . 4376E 00 
0 .6 334E-01 
0 . 5000E-02 
0 .5215E-04 
0.0000 
0 . 1250E-01 
-0 .2500E-02 
0.0000 
0 . 157 IE 01 
0.1571E 01 
0.1323E 08 
0.8684E 06 
0.8684E 06 
0.2166E 07 
0.9537E 06 
0 . 2307E 07 
0 . 456 IE 06 
0.7551E 06 
0.7551E 06 
0 . 16 0 3E-05 
0 . 16 0 IE- 0 4 
0 . 1601E-04 
0.1352E 03 
0 . 23 0 5E 01 
0.2305E 01 
0 . 1 943E 00 
0 . 1144E 08 
0 . 1696E 07 
0 . 1696E 07 
0.4561E 06 
0.7551E 06 
0.7551E 06 
0.2663E 00 
0 . 3947E-01 
0 . 2663E 00 
0 . 3947E-01 


0 . 1000E-01 
0.5580E 00 
0.5524E 00 
0.4420E 00 
0.4376E 00 
0.6334E-01 
0 . 5000E-02 
0 . 5215E-04 
0.0000 
0 .1750E-01 
0 .2500E-02 
0 . 0000 
0 . 157 IE 01 
0.1571E Cl 
0.1323E 08 
0.8684E 06 
0.8684E 06 
0.2166E 07 
0.9537E 06 
0.2307E 07 
0 ♦ 456 IE 06 
0.7551E 06 
0.7551E 06 
0 . 16 0 3E-05 
0 . 1601E-04 
0 . 160 IE-04 
0.1352E 03 
0.2305E 01 
0.2305E 01 
0 . 1943E 00 
Q.1144E 08 
0 . 1696E 07 
0.1696E 07 
0 . 456 IE 06 
0.7551E 06 
0 . 7551E 06 
0.2663E 00 
0 . 3947E-01 
0.2663E 00 
0 . 3947E-01 


0 .2000E-01 
0.5500E 00 
0 . 5390E 00 
0 . 450 0 E 00 
0 . 441 0E 00 
0 . 5443E-01 
0 . 1Q0OE-01 
0 . 5850E-04 
0.0000 
0 .2500E-01 
0 . 1Q0OE-01 
0 .0 0 00 
0.0000 
0.0000 
0.2101E 08 
0.7797E 06 
0.7797E 06 
0.1631E 07 
0.8713E 06 
0.1776E 07 
0 . 3238E 06 
0 . 547 0E 06 
0 . 547 0E 06 
0 . 1418E-06 
0 . 2464E-04 
0 .2464E-04 
0.3195E 03 
0.3702E 01 
0.3702E 01 
0 . 1991E 00 
0 . 1726E 08 
0 . 1127E 07 
0 . 1127E 07 
0.3238E 06 
0 . 547 0 E 06 
Q.5470E 06 
0.2945E 00 
0 . 1922E-01 
0.2945E 00 
0 . 1922E-01 


41 

NUL23 

42 

NUL32 

43 

DPL1 

44 

DPL2 

45 

DPL3 

4 6 

3TAL1 

47 

BTAL2 

48 

ETAL3 

49 

ILMFC 

50 

TEMPD 

51 

LSC11T 

52 

LSC11C 

53 

LSC11D 

54 

LSC22T 

55 

LSC22C 

56 

LSC12 

57 

LSC23 

58 

LSCC23 

59 

LSCC13 

60 

LSCDF 

61 

KL12AB 

62 

MDEIE 

63 

RELP.OT 

64 

EPS11 

65 

EPS22 

66 

EPS12 

67 

SIG11 

68 

SIG22 

69 

SIG12 

70 

DELFI 

71 

HFC 

72 

MPCTGE 

73 

SIG13 

74 

SIG23 

75 

SIG33 


0 . 482 IE 00 
0 . 482 IE 00 
0 .86 0QE-04 
0 . 5168E-04 
0 . 5168E-04 
0 . 4981E-04 
0 . 1452E-02 
0 . 1452E-02 
0.0 0 0 0 
0.0 00 0 
0.2228E 06 
0.8764E 05 
0.8764E 05 
0 . 5006E 04 
0 . 1502E 05 
0.5126E 04 
0 . 3983E 04 
0.0000 
0 .0 0 0 0 
0.0000 
0 . 9858E 00 
0 . 9646E 0 0 
0.0000 
0 .2751E-02 
-0 .2236E-03 
0 . 9946E-09 
0 . 4769E 05 
0 . 6647E 03 
0 . 5441E-03 
0.0000 
0.1121E 01 
0.0000 
0.0000 
0.0000 
0.0000 


0 . 3985E 0 0 
0.3935E 00 
0 .8736E-04 
0 . 5117E-04 
0 .5117E-04 
0 . 9858E-04 
0 . 8565E-03 
0 . 1455E-02 
0.8405E 02 
0.0000 
0 . 2168E 06 
0 . 1665E 06 
0 . 16 65E 06 
0 . 9915E 04 
0 . 2314E 05 
0 . 1195E 05 
0 . 1019E 05 
0 .6196E 05 
0.7303E 05 
0 . 4164E-03 
0 . 9075E 00 
0 . 7801E 00 
0 . 1000E 01 
-0 . 2236E-03 
0 . 2751E-02 
-0 .8536E-08 
-0.1329E 04 
0.4614E 04 
-0 .6445E-02 
-0 . 4765E-08 
0.6393E 00 
0.0000 
0.0000 
0.0000 
0.0000 


0 . 3985E 00 
0.3985E 00 
0 . 8736E-04 
0 .5117E-04 
0 .5117E-04 
0 . 9853E-04 
0 .8565E-03 
0 . 1455E-02 
0.8916E 02 
0.0 0 0 0 
0 . 2168E 06 
0 . 1665E 06 
0 . 1665E 06 
0 . 9915E 04 
Q.2314E 05 
Q.1195E 05 
Q.1019E 05 
0.1351E 06 
0.8147E 05 
0 . 3925E-03 
0 . 9075E 00 
0 . 7801E 00 
0.1000E 01 
-0 . 2236E-03 
0 .2751E-02 
-0 .8536E-08 
-0.1329E 04 
0.4614E 04 
-0 .6445E-02 
0.0000 
0 . 6393E 00 
0.0000 
0.0000 
0.0000 
0 .0 00 0 


0 . 4821E 00 
0.4821E 00 
0 .86 00E-04 
0 . 5168E-04 
0 . 5168E-04 
0 . 4 93 IE- 04 
0 . 1452E-02 
0 . 1452E-02 
0.8405E 02 
0 . 0000 
0.2228E 06 
0.3764E 05 
0.8764E 05 
O.5006E 04 
0 . 1502E 05 
0.5126E 04 
0 . 3983E 04 
0 . 4417E 05 
0.5238E 05 
0 . 4164E-03 
Q.985CE 00 
0 . 9646E 00 
0 . 1000E 01 
0 .2751E-02 
-0 . 2236E-Q 3 
0 . 9946E-09 
0.4769E 05 
0 . 6647E 03 
0 . 5441E-03 
0 . 4765E-08 
0.1121E 01 
0 .0000 
0 .0000 
0 .0000 
0 .0000 
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Item 13 


DERAILS OF POISSON RATIO MISMATCH 


FOISSON"S RATIOS OF THE COMPOSITE 

AHUCXY = 0.0313 

AMUCYX = 0.0303 

AKUCSX = -0.0000 

ANUCSY = 0.0000 


NO. 

THETA 

AMULXY 

ANULSX 

ANULSY 

POIDFN 

POIDFS 

1 

0 . 0 

0.2945 

0.0 0 0 0 

0 . 000 0 

0 .2132 

-0.0 0 0 0 

2 

90 . 0 

0.0395 

-0 . 00 0 0 

-0 . 000 0 

-0 . 0413 

-0 . 0 0 0 0 

3 

90 . 0 

0 .0395 

-0.00 0 0 

-0.00 0 0 

-0 .0416 

-0 . 0 0 00 

4 

0 . 0 

0 . 2945 

0.00 0 0 

0 . 0000 

0.2132 

-0.0 0 00 


Item 14 





FREE 

EDGE STRESSES 



PLY 

THETA 

SIGXX 

SIGYY 

SIGXY YDCAY LENGTH SIGZY 

SIGZZ 

SIGZX 


1 

0 . 0 

0 . 143E 

01 

0 . 199E-01 

0 . 163E-07 

0.119E 00 

0.116E-01 

0 . 337E-02 

0 .412E-08 

1 

0 . 0 

0 . 143E 

01 

0 . 199E-01 

0 . 163E-07 

0 . 000 

0 . 116E-01 

0 . 337E-02 

0 .412E-08 

2 

90 . 0 

0 . 138E 

00 

-0 . 399E-01 

-0 . 326E-07 

0.298E 00 

-0 . 97 9E-08 

0 .803E-03 

0 . 243E-03 

2 

90 . 0 

0 . 138E 

00 

-0 . 399E-01 

-0 . 326E-07 

0 . 345E 00 

-0 .845E-08 

0 . 200E-03 

0.248E-08 

3 

90 . 0 

0 . 138E 

00 

-0 . 399E-01 

-0 . 326E-07 

0.298E 00 

-0 . 574E-08 

0 .803E-03 

0 . 248E-08 

3 

90 . 0 

0 . 1 38E 

00 

-0 . 399E-01 

-0 . 326E-07 

0.345E 00 

-0 . 495E-08 

0 . 200E-03 

0 .248E-08 

4 

0 . 0 

0 . 143E 

01 

0 . 199E-01 

0 . 163E-07 

0.119E 00 

0 . 116E-01 

0 . 337E-02 

0 . 412E-Q8 

4 

0 . 0 

0 . 1 43E 

01 

0 . 199E-01 

0 . 163E-07 

0.000 

0 . 116E-01 

0 . 337E-02 

0 . 412E-08 


MOTE: THE INTERLAMINAR STRESSES ARE BETWEEN PLIES CI-1) AND Cl). 

NOTE: IF THE PLY NO IS REPEATED THEN THE SECOND ONE INDICATES STRESSES IN THE SECONDARY COMPOSITE. 
NOTE: FOR ANGLE PLY LAMINATES SIGYY IS 0. CONSEQUENTLY SIGZY AN D SIGZZ ARE COMPUTED AS ZERO. 

TO OBTAIN NONTRIVIAL SIGZY AND SIGZZ, ONE MUST SPECIFY A THIN INTERPLY LAYER. 

THE INTERPLY LAYER THICKNESS MAY BE OBTAINED FROM THE PLY PROPERTY TABLE. 
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Item 15(a) 

microstresses 


0. o. 
0 . 0 . 
o. 0. 

HOTE^ :""N0 MOISTURE OR TEMPERATURE 
(NOTE: R0I45-PR0PERTY , CO L UMNS-L AYER ) 


FOR LOAD CONDITIONS 
MEMBRANE LOADS NE S ( X , Y , XY-M ) 
BENDING LOADS MBS (X , Y , XY-M) 
av7 ay 7 Awn APPI TFD PRESSURES 


ARE 

ARE 


10 0 0 . 

0. 


MATERIAL^SYSTEM as— /imls 5 £ LA/| j™| 

AS— /IMLS AS — /IMHS 

ORIENTATION 0-0 90 *?_ 


3 

SGLA/HMHS 
AS- - / IMHS 
90 . 0 


A 

AS — /IMLS 
AS--/ IMLS 
0 . 0 


1 

SMIL 

1 

SMIL 

2 

SM1T 

2 

SM1T 

3 

SF1L 

3 

SF1L 

4 

SPIT 

4 

SF1T 

5 

SM2AL 

5 

SM2AL 

6 

SM2AT 

6 

SM2AT 

7 

SM2BL 

7 

SM23L 

8 

SM2BT 

8 

5M23T 

9 

SF2BL 

9 

SF2BL 

10 

SF2BT 

10 

5F2BT 

11 

SM3AL 

11 

SM3AL 

12 

SM3AT 

12 

SM3AT 

13 

SM3BL 

13 

5M3BL 

14 

5M3BT 

14 

SM3BT 

15 

SF3BL 

15 

SF3BL 

16 

SF3BT 

16 

SF3BT 

17 

SM12A 

17 

SM12A 

18 

SM12B 

18 

SM12B 

19 

SF12B 

19 

SF12B 

20 

SM13A 

20 

SM13A 

21 

SM13B 

21 

SM13B 

22 

SF13B 

22 

SFi 33 

23 

SM23A 

23 

SM23A 

24 

SM23B 

24 

SM23B 

25 

SF23B 

25 

SF23B 


0.13S1E 04 
0.0000 
0.2669E 03 
0.0000 
0.8564E 05 
0 . 00 0 0 
-0.2185E 03 
0.0000 
0.1595E 03 
0 . 0000 
0.3445E 03 
0.0000 
-0.1662E 05 
0.0000 
0.7763E 03 
0.0000 
-0.1662E 05 
0.0000 
0.7763E 03 
0.0000 
0.1595E 03 
0.0000 
-0.2126E 02 
0.0000 
-0.1662E 05 
0 .0000 
0 . 23 16 E 03 
0.0 0 0 0 
-0.1662E 05 
0 . 0 0 00 
0.2316E 03 
0 . 0000 
0 .2137E-03 
0 .0000 
0 .6592E-03 
0.0000 
0. 6592 E- 03 
0.0000 
0 .2137E-03 
0.0000 
0 .6592E-03 
0.0000 
0 . 6592E-03 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


-0.8714E 02 
-0.5809E 02 
Q.1534E 04 
0.1561E 04 
-0.1441E 04 
-0 . 36 0 2E 04 
-0.4087E 03 
- 0 . 24 0 6 E 04 
-0.7294E 01 
-0.4862E 01 
0 . 17 0 6 E 04 
0 . 2 324 E 04 
0.4632E 03 
0.4314E 03 
0.3858E 04 
0.7882E 04 
0 .4632E 03 
0.4314E 03 
0.3858E 04 
0.7882E 04 
-0.7294E 01 
-0 . 486 2 E 01 
-0.9887E 02 
-0.6591E 02 
0 . 46 32 E 03 
0.4314E 03 
0. 1607E 04 
0 . 1497E 04 
0.4632E 03 
0.4314E 03 
0.1607E 04 
0.1497E 04 
-0 .2348E-02 
-0 .2439E-Q2 
-0 .6382E-C2 
-0 . 9945E-02 
-0.6382E-02 
-0 . 9945E-02 
- 0 . 2348 E- 0 2 
-0.2439E-02 
-0 .6382E-02 
-0 . 9945E-02 
-0 . 6382E-02 
-0 . 9945E-02 
0.0000 
0.0000 
0.0000 
0 .00 0 0 
0.0000 
0.0000 


-0 .8714E 02 
-0.5809E 02 
0.1534E 04 
0. 1561E 04 
-0.1441E 04 
-0 .3602E 04 
-0.4087E 03 
-Q.2406E 04 
-0.7294E 01 
-0.4862E 01 
0.1706E 04 
0 .2324E 04 
0.4632E 03 
0.4314E 03 
0.3858E 04 
0.7882E 04 
0.4632E 03 
0 .4314E 03 
0.3358E 04 
0 . 7882 E 04 
-0.7294E 01 
-0.4862E 01 
-0.9887E 02 
-0 . 6 5 9 1 E 02 
0 . 46 32 E 03 
0.4314E 03 
0.1607E 04 
0 . 1497E 04 
0.4632E 03 
0 .4314E 03 
0 . 16 0 7 E 04 
0.1497E 04 
-0 . 2 348 E- 02 
-0 .2439E-02 
-0 .6382E-02 
-0.9945E-02 
-0.6 382 E- 0 2 
-0 . 9945E-02 
-0 . 2 348 E- 02 
-0.2439E-02 
-0 . 6382E-02 
-0 . 9945E-02 
-O . 6382E-02 
-0 . 9945E-02 
0.0000 
0.0000 
0.0000 
0 .0000 
0 . 0000 
0 . 000 0 


0 . 1381E 04 
0 . 0000 
0.2669E 03 
0.000 0 
0.8564E 05 
0.0000 
-0.2185E 03 
0 . 0000 
0.1595E 03 
0 .00 00 
0.3445E 03 
0 . 0000 
-0.1662E 05 
0.0000 
0.7763E 03 
0.0000 
-0.1662E 05 
0.0000 
0.7763E 03 
0.0000 
0.1595E 03 
0.000 0 
-0.2126E 02 
0.0000 
-0.1662E 05 
0 . 0000 
0.2316E 03 
0 .0000 
-0.1662E 05 
0 . 0000 
0.2316E 03 
0 .0000 
0.2137E-03 
0 .0000 
0 . 6592E-03 
0 .0000 
0.6592E-O3 
0 .0000 
0 .2137E-03 
0 .0000 
0 .6592E-03 
0 .0000 
O . 6592E-03 
0.0000 
0.0000 
0.0000 
0 .0 000 
0 .0000 
0 . 000 0 
0.0000 


NOTATION: S STRESS (SIGMA) f . Tneo 

M — - MATRIX AND F — FIBER ayc* 

l 2 ,3 — DIRECTIONS FOR STRESSES - PLY MATERIAL AXES 

LfT DIRECTIONS OF PLY STRESSES 

A — - REGION CONTAINING NO FIBERS 
B REGION CONTAINING FIBERS AND MATRIX 

FYAMPLF: SM2AL STANDS FOR TRANSVERSE NORMAL STRESS 

EXAMPLE- T0 A L0AD i N THE LOGITUDINAL 

DIRECTION 
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Item 15(b) 

MICROSTRESS INFLUENCE COEFFICIENTS 


THE FOLLOWING ARE THE MICROSTRESS INFLUENCE COEFFICIENTS FOR THE PRIMARY COMPOSITE AS— /IMLS SYSTEM 


INF. COEF. 

SIGMA 11 
LSS/SQ. IN 

SIGMA22 
LtiS/SO. IN 

5IGMA12 
LBS/SO. IN 

SIGMA 1 3 
LBS/SQ. IN 

SIGMA 2 3 
LBS/SQ. IN 

DELTA T 
1 DEG F 

DELTA M 
IX 

1 SMI I 

0.0290 

0 .AC15 

0 . 0000 

0.0000 

0.0000 

-28.A291 

-1975.0911 

2 SM22A 

0 .0033 

0.5133 

0.0000 

0.0000 

0.0000 

-16 . 1 8 1 A 

-127A. 0061 

3 SM22B 

-0 . 3 ASA 

1 . 1678 

O.OCOO 

0.0000 

0.0000 

5.6376 

AA3.8633 

A SMI 2 A 

0.0000 

0 . OCOO 

0.3927 

0 . 0000 

0.0000 

0.0000 

0.0 0 0 0 

5 SM12B 

0 .0 0 0 0 

0 . 0 0 0 0 

1.2116 

0 . 0000 

0 . 0000 

0.0000 

0.0 0 0 0 

6 SM13A 

0 . 0 0 0 0 

0 . 00 00 

0 . 0000 

0.3927 

0 . 0000 

0.0000 

0.0000 

7 SM13B 

0 . 0 0 0 0 

0 . OCOO 

0 . 0000 

1.2116 

0.0000 

0.0000 

0.0 000 

8 SM23A 

0.0000 

0 . OCOO 

0 .0000 

0.0000 

0.5A75 

0 . 0 0 0 0 

0.0 000 

9 SM23B 

0 .0000 

0 . OCOO 

0 . 0000 

0.0000 

1 .A0A3 

0 .0000 

0 . 0000 

10 SM3 3 A 

0 . 0033 

-0 . 0320 

0 .0000 

0.0000 

0.0000 

-16 . 181 A 

-127 A .0 06 1 

11 SN33B 

-0 . 3 AS A 

0 . 3 A 8 A 

0 . 0000 

0.0000 

0 .0000 

5.6376 

AA3.8633 

12 5 FI 1 

1.7955 

-0.3258 

0 . 0000 

0.0000 

0.0000 

21 . A A 6 5 

1 5 AA .1631 

13 SF22B 

-0. 3ASA 

1 . 1678 

0 .0000 

0.0000 

0.0000 

5.6376 

AA3 .8633 

14 SF33B 

-0 . 3 ASA 

0 . 3ASA 

0 .0000 

0.0000 

0.0000 

5.6376 

AA3 .8633 

IS SF12 

0.0 00 0 

0.0000 

1.2116 

0.0000 

0.0000 

0 . 0000 

0 . 0000 

16 S FI 3 

0.0000 

0 .0000 

0 .0000 

1.2116 

0.0000 

0 . 0000 

0 . 0 0 0 0 

17 SF23B 

0.000 0 

0 .0000 

0.0000 

0.0000 

1.AOA3 

0 .0000 

0 .0 000 

NOTE: TO OBTAIN THE ABSOLUTE VALUE OF THE 

: MICROSTRESSES THE INF. COEF. SHOULD BE MULTIPLIED BY THE APPROPRIATE 

STRESSES OR THE TEMPERATURE 

GRADIENT OR THE MOISTURE CONTENT. 




EXPLANATION: SM22B 

.FOR EXAMPLE, 

STANDS FOR TRANSVERSE NORMAL STRESS INFLUENCE 

COEFFICIENT IN REGION B 


M 

I C R 0 S T R 

ESS I N F L U 

E N C E COE 

F F I C I E N 

T S 


THE FOLLOWING ARE THE MICROSTRESS INFLUENCE COEFFICIENTS 

FOR THE PRIMARY 

COMPOSITE SGLA/HMHS SYSTEM 

INF. COEF. 

5IGMA11 

5IGMA22 

SIGMA 1 2 

SIGMA13 

SIGMA23 

DELTA T 

DELTA M 


LBS/SQ . IN 

LBS/SQ . IN 

LBS/SO. IN 

LBS/SQ. IN 

LBS/SQ. IN 

1 DEG F 

IX 

1 SMI I 

0 . 0655 

0.3325 

0.0000 

0.0000 

0.0000 

-28.7981 

-2926 .0596 

2 Si iZ 2 A 

0.0055 

0.3655 

0.0000 

0.0000 

0.0000 

-17 . 9927 

-2357 .6A9A 

3 5M22B 

-0 . 3AS A 

0 .S363 

0 .0000 

0.0000 

0.0000 

6 .2687 

821 . A0A3 

A 5M12A 

0.0 0 0 0 

0.0000 

0 . 36 A 3 

0 . 0000 

0 . 0000 

0.0 0 00 

0.0000 

5 SM12B 

0 . 0000 

0 . 0 0 0 0 

0 .9902 

0.00 00 

0.0000 

0.0000 

0.0000 

6 5M13A 

0.0000 

0 .0 0 0 0 

0 .0000 

0 . 36 A3 

0.0000 

0.0000 

0.0000 

7 SMI 3 3 

0. 0000 

0.0000 

0.0000 

0. 9902 

0.0000 

0 .0000 


3 S 1 1 Z 3 A 

0.0 0 0 0 

0.0000 

0.0000 

0 . 0 000 

0.6091 

0.0000 

0.0000 

9 SM233 

0.000 0 

0.0000 

0 .0000 

0.0000 

2.0A23 

0.0000 

0.0000 

10 SM33A 

0 . 0 055 

-Q . 021A 

0.0000 

0.0000 

0.0000 

-1 7 . 9927 

-2357 . 6A9A 

11 ^ M 3 3 3 

-0 . 3 AS A 

0 . 3A8A 

0 .0000 

0.0000 

0.0000 

6.2687 

821 . AQA3 

12 S r 1 1 

1 . 0S37 

-0 . 0386 

0 .0000 

0 .0000 

0.0000 

-1A.8A8A 

1222 .AA17 

13 3F22B 

-0 . 3 ASA 

0,5363 

0 . 0 0 0 0 

0 . 0000 

0.0000 

6.2687 

821 . A 0 A 3 

i -* j F 3 3 u 

- 0 . 3 AS A 

0 . 3 AS A 

0 .0000 

0.0000 

0.0000 

6.2687 

821 . A0A3 

15 5 FI 2 

0.0 0 00 

0.0000 

0 .9902 

0.0000 

0.0000 

0.0000 

0 . 0 0 0 0 

16 5 F 1 3 

0.0000 

0 . 0000 

0 . 0000 

0.9902 

0.0000 

0 .0000 

0 . 0000 

17 3F2 3B 

0.0000 

0 . 00 0 0 

0 . 0000 

0.0000 

2. 0A23 

0.0000 

0.0000 

NOTE: TO OBTAIN THE ABSOLUTE VALUE OF THE 

MICROSTRESSES THE INF. COEF. SHOULD BE MULTIPLIED BY THE APPROPRT ATF 

SIKcSSES OR THE TEMPERATURE 

GRADIENT OR THE MOISTURE CONTENT. 




EXPLANATION: SM22B- 

, FOR EXAMPLE, 

STANDS FOR TRANSVERSE NORMAL STRESS INFLUENCE 

COEFFICIENT IN 

REGION B 




MICROSTRE 

S 5 I N F L U 

ENCE COEF 

F I C I E 

NTS 


THE 

FOLLOWING 

ARE THE MICROSTRESS INFLUENCE 

COEFFICIENTS 

FOR THE SECONDARY 

COMPOSITE 

AS--/IMHS SYSTEM 


INF. 

COEF. 

SIGMA11 SIGMA22 

SIGMA 12 

SIGMA 13 

SIGMA 2 3 

DELTA T 

DELTA M 




LBS/SQ. IN LBS/SQ. IN 

LBS/SQ. IN 

LBS/SQ. IN 

LBS/SQ. IN 

1 DEG F 

IX 


1 SMI 1 

0 .QA37 

0 .338A 

0 . oooo 

0.0000 

0.0000 

-17 .1987 

-1 950 .7063 

2 SM22A 

0.0037 

0.5036 

0.0000 

0.0000 

0.0000 

-9.9952 

-157 1 . 766 1 

3 SM22B 

-0 . 32A5 

1.7084 

0 . 0000 

0.0000 

0.0000 

3.2A38 

510 . 0889 

A S<T12A 

0.0 0 0 0 

0.0000 

0. 378A 

0 .oooo 

0.0000 

0 .0000 

0.0000 

3 SM12B 

0.0000 

0.0000 

1 .5A30 

0 .0000 

0 .0000 

0.0000 

0 . oooo 

6 SM13A 

0 . 0 0 0 0 

0 . 0000 

0 . 378A 

0.0000 

0.0000 

0 .0000 

0 . oooo 

7 S III 33 

0.0 0 0 0 

0.0000 

1 . 5 A 3 0 

0.0000 

0 .0000 

a . oooo 

0.0000 

8 Si 12 3 A 

0 . 0000 

0.0000 

0 . OOOO 

0.0000 

0.4060 

0.0000 

0.0000 

9 SM23B 

0.0000 

0 .0000 

0 .0000 

0 .0000 

1.0551 

0 .0000 

0.0000 

10 SM33A 

0.0037 

-0 . 01A3 

0 .0000 

0.0000 

0.0000 

-9 . 9952 

-1571 .7661 

11 SM33B 

-0 . 32A5 

0.32A5 

0.0000 

0 .0000 

0.0000 

3. 2A38 

510 . 0389 

12 SF1 1 

2.70 93 

-0.5215 

0.0000 

0.0000 

0.0000 

66 . 7289 

3056 .1045 

1 3 SF22B 

-0 . 32 A 5 

1 .7 08A 

0.0000 

0.0000 

0.0000 

3. 2A38 

510 . 0889 

1A SF33B 

- 0 . 32A5 

0.32A5 

0 . 0000 

0 . 0000 

0.0000 

3. 2A38 

510 . 0889 

15 5F12 

0.0 0 0 0 

0 .0 00 0 

1 .5430 

0.0000 

0.0000 

0.0000 

0.0000 

16 SF1 3 

0.0000 

0.0 00 0 

1 . 5 A 3 0 

0.0000 

0.0000 

0.0000 

0 . oooo 

17 SF23B 

0.0 0 0 0 

0 . 00 00 

0 . 0 0 0 0 

0.0000 

1.0551 

0 .0000 

0.0000 


NOTE: TQ OBTAIN THE ABSOLUTE VALUE OF THE MICROSTRESSES THE INF. COEF 

STRESSES OR THE TEMPERATURE GRADIENT OR THE MOISTURE CONTENT. 


SHOULD BE MULTIPLIED BY THE APPROPRIATE 


EXPLANATION'. SM22B , FOR EXAMPLE, STANDS FOR TRANSVERSE NORMAL STRESS INFLUENCE COEFFICIENT IN REGION B 
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NOTE : K1XX — > STRESS CONCENTRATION FACTOR DUE TO SIGMA XX 

K1YY — > STRESS CONCENTRATION FACTOR DUE TO SIGMA YY 

K1XY — > STRESS CONCENTRATION FACTOR DUE TO SIGMA XY 

LAYUP --> 0 90 90 0 


THETA K1XX K1YY K1XY THETA K1XX K1YY KIXY 


0.0 

-0.6160 

5.0 

-0.5709 

10.0 

-0.4572 

15.0 

-0.3168 

20.0 

-0.1799 

25.0 

-0 . 0569 

30.0 

0.0532 

35.0 

0 . 1566 

40 . 0 

0.2613 

45 . 0 

0.3764 

50.0 

0.5138 

55.0 

0.6398 

60 . 0 

0 . 9302 

65.0 

1.2764 

70 . 0 

1 .7980 

75.0 

2.6025 

80 . 0 

3.7983 

85.0 

5.2286 

90 . 0 

5.9765 

95.0 

5.2287 

100.0 

3.7984 

105.0 

2.6026 

110.0 

1.7980 

115.0 

1.2764 

120.0 

0 . 9302 

125.0 

0 . 6899 

130.0 

0 .5138 

135.0 

0.3764 

140 . 0 

0.2613 

145.0 

0.1567 

150.0 

0.0532 

155.0 

-0.0569 

160 . 0 

-0 . 1799 

165.0 

-O . 3168 

170.0 

-0 . 457 1 

175.0 

-0 . 57 0 9 


3.8562 0.0000 

3.6729 -1.1975 

3.2156 -2.1209 

2.6650 -2.6886 

2.1516 -2.9777 

1.7253 -3.1020 

1.3875 -3.1493 

1.1225 -3 . 1741 

0.9117 -3.2080 

0.7383 -3.2701 

0.5879 -3.3730 

0.4473 -3.5263 

0.3019 -3.7367 

0.1331 -A. 0029 

-0.0864 -4.2958 

-0.3968 -4.5006 

-0.8378 -4.2743 

-1.3548 -2.9009 

-1.6233 -0.0001 

-1.3549 2.9007 

-0.8379 4.2743 

-0.3969 4.5007 

-0.0865 4.2958 

0.1330 4.0029 

0.3018 3.7367 

0.4472 3.5263 

0.5879 3.3730 

0.7382 3.2702 

0.9116 3.2081 

1.1225 3.1741 

1.3875 3.1493 

1.7252 3.1021 

2.1515 2.9777 

2.6648 2.6886 

3.2155 2.1210 

3.6728 1.1976 


180.0 

-0.6160 

185.0 

-0.5709 

190 . 0 

-0.4572 

195.0 

-0.3168 

200 . 0 

-0.1799 

205.0 

-0.0570 

210.0 

0.0532 

215.0 

0.1566 

220 . 0 

0.2613 

225.0 

0.3764 

230.0 

0.5137 

235.0 

0.6898 

240.0 

0.9301 

245.0 

1.2763 

250.0 

1.7979 

255.0 

2.6023 

260.0 

3.7980 

265.0 

5.2283 

270.0 

5.9765 

275.0 

5.2290 

280 . 0 

3.7987 

285.0 

2.6028 

290.0 

1.7982 

295.0 

1.2765 

300.0 

0.9303 

305.0 

0.6899 

310.0 

0.5138 

315.0 

0.3764 

320.0 

0.2613 

325.0 

0.1567 

330.0 

0.0532 

335.0 

-0.0569 

340.0 

-0.1798 

345 . 0 

-0.3168 

350 . 0 

-0.4571 

355.0 

-0.5709 


3.8562 0.0002 

3.6729 -1.1973 

3.2155 -2.1208 

2.6649 -2.6885 

2.1515 -2.9777 

1.7252 -3.1020 

1.3875 -3.1493 

1.1225 -3.1741 

0.9116 -3.2080 

0.7382 -3.2701 

0.5879 -3.3730 

0.4472 -3.5262 

0.3018 -3.7367 

0.1331 -4.0028 

-0.0865 -4.2957 

-0.3969 -4.5006 

-0.8379 -4.2744 

-1.3549 -2.9013 

-1.6233 -0.0008 

-1.3550 2.9003 

-0.8380 4.2741 

-0.3970 4.5007 

-0.0865 4.2958 

0.1330 4.0030 

0.3018 3.7368 

0.4472 3.5263 

0.5878 3.3731 

0.7382 3.2702 

0.9116 3.2081 

1.1224 3.1741 

1.3874 3.1493 

1.7251 3.1021 

2.1514 2.9778 

2.6647 2.6887 

3.2153 2.1211 

3.6728 1.1979 
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Item 17 

LOCATIONS OF PROBABLE DELAMINATION 


RESULTS FOR PLY NO. 1 ORIENTATION 0 

.0 MATERIAL AS- 

-IMLS AS — IMLS 



CRITERION 

RANGE 

VALUE 

LOCATION x 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

0.0— 90.0 

0.588 

75 0 

MAX 

OF 

KIYYM(NUCRT-NULRT) 

0.0-- 90.0 

0 .822 

0 0 

MAX 

OF 

KIXYx(NUCRT-NULRT) 

0.0— 90.0 

1 . All 

60.0 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

90.0— 180.0 

0 . 588 

105.0 

MAX 

OF 

K1YYMNUCRT-NULRT) 

90.0— 180.0 

0.822 

180 . 0 

MAX 

OF 

K1XYMNUCRT-MULRT) 

90.0— 180.0 

1 . All 

120 . 0 

MAX 

OF 

K IXX# ( NUCRT-NULRT) 

180.0— 270.0 

0 . 588 

255 . 0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

180. 0-- 270.0 

0.822 

180.0 

MAX 

OF 

K1XY#( NUCRT-NULRT ) 

180.0— 270.0 

1 .All 

2A0 . 0 

MAX 

OF 

K1 XX ^(NUCRT-NULRT) 

270.0— 0 .0 

0.588 

285 . 0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

270.0— 0.0 

0.822 

0 0 

MAX 

OF 

K1XY*( NUCRT-NULRT) 

270. 0-- 0.0 

l.AU 

300 . 0 

RESULTS FOR PLY NO. 2 ORIENTATION 90. 

.0 MATERIAL SGLAHMHS AS — IMHS 



CRITERION 

RANGE 

VALUE 

LOCATION x 

MAX 

OF 

K1XX#<NUCRT-NULRT) 

0.0— 90.0 

1 . A07 

90.0 

MAX 

OF 

KIYYK(NUCRT-NULRT) 

0.0-- 90.0 

1 . 087 

15 0 

MAX 

OF 

K1XY*(NUCRT-NULRT) 

0.0-- 90.0 

I.A76 

25 . 0 

MAX 

OF 

K1XX#(NUCRT-NULRT) 

90.0 — 180.0 

1.A07 

90 . 0 

MAX 

OF 

K1YY*(NUCRT-NULRT) 

90. 0— 180.0 

1.087 

165 . 0 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

90. 0-- 180.0 

1 . A76 

155.0 

MAX 

OF 

KlXXx (NUCRT-NULRT) 

180. 0-- 270.0 

1 . A07 

270.0 

MAX 

OF 

KIYY#<NUCRT-NULRT) 

180.0— 270.0 

1 . 087 

195.0 

MAX 

OF 

K1XY#<NUCRT-NULRT) 

180.0— 270.0 

1 . A76 

205 . 0 

MAX 

OF 

K IXX* ( NUCRT-NULRT) 

270.0— ♦ 0.0 

1. A07 

270.0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

270. 0-- 0.0 

1.087 

3A5. 0 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

270. 0-- 0.0 

1.A76 

335.0 

RESULTS FOR PLY NO. 3 ORIENTATION 90. 

0 MATERIAL SGLAHMHS AS — IMHS 



CRITERION 

RANGE 

VALUE 

LOCATION x 

MAX 

OF 

KiXX#(NUCRT-NULRT) 

0.0— 90.0 

1 . A07 

90 . 0 

MAX 

OF 

KIYYM(NUCRT-NULRT) 

0.0— 90.0 

1 . 087 

15 . 0 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

0.0— 90.0 

1 . A76 

25 . 0 

MAX 

OF 

K IXX* (NUCRT-NULRT ) 

90 . 0— 180.0 

1 . A07 

90 . 0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

90.0— 180.0 

1 . 087 

165.0 

MAX 

OF 

KIXYXCNUCRT-NULRT) 

90. 0-- 180.0 

1 . A76 

155.0 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

180. 0-- 270.0 

1.A07 

270 . 0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

180.0— 270.0 

1.087 

195. 0 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

180.0— 270.0 

1. A76 

205 . 0 

MAX 

OF 

KIXXM(NUCRT-NULRT) 

270.0— 0 .0 

1. A07 

270.0 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

270.0— 0.0 

1.087 

3A5 . 0 

MAX 

OF 

KIXYM(NUCRT-NULRT) 

270.0— 0.0 

1.A76 

335.0 

RESULTS 

FOR PLY NO. A ORIENTATION 0 . 

0 MATERIAL AS— 

IMLS AS — IMLS 



CRITERION 

RANGE 

VALUE 

LOCATION x 


MAX 

OF 

KIXXX(NUCRT-NULRT) 

0.0 — 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

0 . o — 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

0 . 0 — 

MAX 

OF 

KIXXX(HUCRT-NULRT) 

90.0 — 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

90,0 — 

MAX 

OF 

K1XYXCNUCRT-NULRT) 

90.0 — 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

180.0 — 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

180.0 — 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

180 . 0 — 

MAX 

OF 

KIXXX(NUCRT-NULRT) 

270 . 0 — 

MAX 

OF 

KIYYX(NUCRT-NULRT) 

270 . 0 — 

MAX 

OF 

KIXYX(NUCRT-NULRT) 

270. Q — 


90.0 

0.588 

75.0 

90.0 

0.822 

0 . 0 

90.0 

1.411 

60 . 0 

180.0 

0.588 

105.0 

180.0 

0.822 

180.0 

180.0 

1.411 

120.0 

270.0 

0.588 

255.0 

270.0 

0.822 

180 . 0 

270.0 

1.411 

2 A 0 . 0 

0 . 0 

0.538 

285.0 

0 . 0 

0.822 

0.0 

0.0 

1. All 

300 . 0 


NOTES ; K1XX — > STRESS CONCENTRATION FACTOR DUE TO SIGMA XX 

K1YY — > STRESS CONCENTRATION FACTOR DUE TO SIGMA YY 

K1XY — > STRESS CONCENTRATION FACTOR DUE TO SIGMA XY 

NULRT — > PLY POISSON RATIO IN R AND T AXES 

NUCRT — > COMPOSITE POISSON RATIO IN R AND T AXES 
( R AND T ARE THE RADIAL AND THE TANGENTIAL DIRECTIONS) 
ONLY 5 DEG. INTERVALS ARE CONSIDERED. THE ACTUAL VALUE 
IS EXPECTED TO BE WITHIN 5 DEG. OF THE PRINTED RESULT. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Item 18 


SYSTEM 
AS — /IMIS 


AS--/IMHS 


3 5GLA/HMH5 
AS — / IP1HS 


90. 


AS--/ IMIS 


STRESS 

AND ST 

RAIN IN 

F L U E N 

C E C O 

E F F I C I 

ENTS A 

R R A Y S 


A RESPONSE 

NX 

(UNIT 

NY 

LOAD ... LB. 

NXY 

/INCH) 

MX MY 

(UNIT MOMENT ... LB. 

MXY 

IN/INCH) 

DELTAT 
(1 DEG F) 

DEITAM 

Cl 

0 

EPS11 

EPS22 

EPS12 

SIGll 

SIG22 

SIG12 

0 

2.7511 
-0 .2236 
0.0000 
47 .6932 
0 .6647 
0.0000 

-0.2236 
7 . 2A 90 
-0 . 0000 
-1. A628 
8.1391 
-0 . 0 0 0 0 

0.0000 
-0 .0 0 0 0 
5A . 0302 
0.0000 
-0 . 0 0 0 1 
29.5831 

266 . 00 17 
-59.2160 
0 . QOOO 
A598 .8906 
21 . 666 A 
0.0000 

-59.2159 
2938.3235 
-0 .0 024 
-47.6819 
3309.6074 
-0.0013 

0 .0000 
-0 . 0024 
8011.9414 
-0.0001 
-0 .0027 
4382.7227 

1 . 1801 
5.8164 
-0.0000 
11.7478 
-20.9790 - 
- 0.0 0 0 0 

105.0192 
342 .4932 
-0 .0023 
588 .3518 
1238 .7244 
-0.0012 

EPS11 
EPS22 
EPS 1 2 
SIGll 
SIG22 
SIG12 

0 

-0.2236 
2.7511 
-0.0000 
-1 .3294 
A. 6136 
-0 .0 0 0 0 

7.2A90 

-0.2236 

0.0001 

83.7219 

2.9257 

0.0001 

0.0000 
-0 .000 1 
-5 A . 0802 
0.0002 
-0 .0001 
-A0 .8337 

-14 .80A0 
66.5005 
-0 .0 002 
-140 . 8320 
107.2302 
-0.0002 

734.5815 
-14.8040 
0.0025 
8487 .6133 
309. 9329 
0.0019 

0.0019 
-0 .0 025 
-2002 . 9873 
0.0213 
-0.0034 
-1512.3735 

5.8164 
1.1801 
0.0 0 0 1 
41 . 9581 
-23.4957 - 
0.0000 

342.4929 
105.0192 
0.0029 
2477 . 4478 
•1176 .7 065 
0 .0022 

EPS 1 1 
EPS22 
EPS12 
SIGll 
SIG22 
SIG12 

-0 . 2236 
2.7511 
-0.0000 
-1.3294 
A .6136 
-0.0000 

7.2490 

-0.2236 

0.0001 

83.7219 

2.9257 

0.0001 

0 .0000 
-0.0001 
-5 A .08 0 2 
0.0002 
-0.0001 
-A0 .8337 

14.8040 
-66 .50 04 
0.00 02 
140.8316 
-107.2300 
0 . 0 002 

-734.5796 
14.8040 
-0.0025 
-8487 .5898 
-309 . 9319 
-0.0019 

-0.0019 
0 .0025 
2002 . 9846 
-0 . 0213 
0 .0034 
1512.3713 

5.8164 
1.1801 
0.0001 
41 . 9580 
-23.4957 - 
0.0000 

342 .4927 
105.0192 
0 .0029 
2477 . 4451 
-1176 . 7065 
0.0022 

0 

EPS 1 1 
EPS22 
EPS 1 2 
SIGll 
SIG22 
SIG12 

2 . 7511 
-0.2236 
0.0000 
A 7 .6 932 
0 . 66A7 
0.0 000 

-0.2236 
7.2A90 
-0 . 0000 
-1 . A628 
8.1391 
-0 . 0000 

0.0 0 0 0 
- 0.0000 
5 A . 0802 
0.00 00 
-0 . 0 0 0 1 
29.5831 

-266 .0017 
59.2159 
-0.0 0 00 
-4598.8906 
-21.6664 
-0 . 0 0 0 0 

59.2159 

-2938.3228 

0.0024 

47.6819 

-3309.6067 

0.0013 

-0.0000 
0.0024 
-8011 . 9414 
0.00 01 
0 .0027 
-4382 . 7227 

1.1801 
5.8164 
- 0.0000 
11.7478 
-20 . 9790 
- 0 . 0 00 0 

105.0192 
342 . 4924 
-0 . 0023 
588 . 3518 
-1238 .7251 
-0.0012 


NOTE : STRAINS ARE IN MICRO INCH/INCH. 
STRESSES ARE IN POUNDS/INCH $Q. 


EXPLANATION OF THE INFLUENCE COEFFICIENTS 


NX, NY AND NXY ARE UNIT LOADS IN LB/INCH. MX, MY AND MXY 
ARE UNIT MOMENTS IN LB. IN/INCH. D^LTAT IS A UNIT TEMP 
DIFF AND DEITAM IS A UNIT PERCENTAGE OF MOISTURE CONTENT. 
TO OBTAIN RESPONSE R FOR A GENERAL APPLIED LOAD 
VECTOR F USE THE FOLLOWING EQUATION: 


(R) = ( AINF ) X (F) 


NOTE: 


R IS A 6X1 COLUMN VECTOR DEFINED BY 

(R) = ( EPS II EPS22 EPS 12 SIGH 5IG12 SIG12) 
F IS A 8X1 COLUMN VECTOR DEFINED BY 


(F) = (NX NY NXY MX MY MXY DELTAT DEL .. _..- c 
AINF IS A (6X8) MATRIX CONTAINING THE INFLUENCE 
COEFFICIENTS ARRAYS . 


PLY STRESS INFLUENCE COEFFICIENTS ARRAYS 


PLY NO. 

MATERIAL 

SYSTEM 

THETA 

RESPONSE 

NX 

(UNIT 

SCALE 

NY 

LOAD ... LB 
FACTOR = 
(1/TC) 

NXY 

./INCH) 

33.333 

MX MY MXY 

(UNIT MOMENT ... LB. IN/INCH) 
SCALE FACTOR = 6666.668 

(6/TC##2> 

DELTAT 
(1 DEG F) 
13.354 

DELTAM 
(1 X) 
1223.138 

1 

AS--/IMLS 

0 . 0 

SIGll 

SIG22 

SIG12 

1.4308 

0.0199 

0.0000 

-0 . 0439 
0.2442 
-0 . 0000 

0.0000 
-0 .0 0 0 0 
0.8875 

0 .6898 
0 .0032 
0.0000 

-0.0072 
0 .4964 
-0.0000 

-0 .0000 
-0 .000 0 
0.6574 

0 .8797 
-1.5709 
-0 . 0000 

0.4810 
-1.0127 
-0 .0000 

2 

SGLA/HMH5 

AS--/IMHS 

90 . 0 

SIGll 

5IG22 

SIG12 

-0.0399 
0.1384 
-0 . 0000 

2.5117 
0 . 0878 
0.0000 

0 .0000 
-0.0000 
-1.2250 

-0 . 0211 
0.0161 
-0 .0000 

1.2731 
0 . 0465 
0.0000 

3.0000 
-0 . 00 0 0 
-0.2269 

3.1419 

-1.7594 

0.0000 

2.0255 
-0 . 9620 
0.0 0 0 0 

3 

SGLA/HMH5 
AS — /IMHS 

90 . 0 

SIGll 

SIG22 

SIG12 

-0 . 0399 
0 . 1384 
-0 .00 00 

2.5117 
0 . 0878 
0.0000 

0.0000 
-0.0000 
-1 .2250 

0.0211 
-0.0161 
0 . 0 0 0 0 

-1.2731 
-0.0465 
-0 . 000 0 

-0.0000 

0.0000 

0.2269 

3.1419 

-1.7594 

0.0000 

2.0255 
-0 . 9620 
0.0000 

A 

AS— /IMLS 

0. 0 

SIGll 

SIG22 

SIG12 

1.4308 
0.0199 
0 .0000 

-0 . 0439 
0.2442 
-0 . 0 0 0 0 

0 .0000 
-0 . 0 0 00 
0.8875 

-0 . 6898 
-0 .0032 
-0 . 0000 

0 .0072 
-0 .4964 
0 .0000 

0 .0000 
0.0 0 0 0 
-0.6574 

0 .8797 0.4810 

-1.5709 -1.0127 

-0 . oooo -o . oooo 


NOTE: 


HE MEMBRANE STRESSES ARE NORMALIZED W.R.T THE 
ERAGE STRESS DUE TO UNIT LOAD IN AN EQUIVALENT 
MOGENEOUS SECTION. THE BENDING STRESSES ARE 
RMALIZED U.R.T THE MAXIMUM STRESS DUE TO UNIT 
iMENT. THE TEMPERATURE AND MOISTURE STRESSES ARE 
IRMAL IZED W.R.T THE AVERAGE STRESSES DUE TO UNIT 
MPER ATURE DIFFERENCE AND UNIT PERCENTAGE OF MOIS- 
IRE. TO OBTAIN THE ABSOLUTE VALUES OF THE 5TR- 
iS ES THE INFLUENCE COEFFICIENTS SHOULD BE MULT I- 
. I ED BY THE INDICATED SCALE FACTORS. THESE SHOULD 
: MULTIPLIED BY THE CORRESPONDING LOADS TO OBTAIN 
: RESSES IN THE PLIES. 
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Item 19 

laminate failure stress analysis 


LAYUP — > o 90 90 0 

STRESSES _BASED^UPCH^FIPoT^ PLY FAILURE CRITERIA (NO TEMPERATURE OR MOISTURE 


STRESSES) 


PLY NO. 

- 1 

THETA = 

0.00 MATERIAL 

SYSTEM = AS— IMLS 

AS— IMLS 



LOADS 

SLUT 

222.7741 

KSI 

SL11C 
87 .6392 
KSI 

SL22T 
5 .0 065 
KSI 

SL22C 

15.0194 

KSI 

SL12S 

5.1261 

KSI 

FAIL. LOAD 
KSI 

MODE 

SCXXT MIN ( 
SCXXC MIN ( 
SCYYT MIN ( 
SCYYC MIN ( 
SCXYS MIN ( 

155.699 
-155.699 
-5076 .305 
5076 . 305 
0 . 000 

-61.252 
61.252 
1997 .015 
-1997 .015 
0 . 000 

251 . 063 
-251 . 063 
20.504 
-20.504 
********** 

-753.188 

753.188 

-61.511 

61.511 

********** 

0.000 ) 
0.000 ) 
0.000 ) 
0.000) 
5.776) 

155.699 

61.252 

20.504 

61.511 

5.776 

SL11T 

SL11C 

SL22T 

SL22C 

SL12S 


FAILURE STRESSES BASED UPON FIRST PLY FAILURE CRITERIA 


(NO TEMPERATURE OR MOISTURE STRESSES) 


PLY NO. 

i 

i 

i 

i 

rj i 
i 
i 
i 

H 1 
1 
1 

THETA = 90. 

.00 MATERIAL 

SYSTEM = SGLAHMHS AS— IMHS 



LOADS 

SLUT 

216.8321 

KSI 

SL11C 

166.5112 

KSI 

SL22T 

9.9151 

KSI 

SL22C 

23.1353 

KSI 

SL12S 

11.9513 

KSI 

FAIL. LOAD 
KSI 

MODE 

SCXXT MIN ( -5436.801 
SCXXC MIN ( 5436.801 
SCYYT MIN ( 56.330 
SCYYC MIN ( -86.330 
SCXYS MIN (********** 

4175.063 

-4175.063 

-66.295 

66.295 

********** 

71.638 

-71.633 

112.967 

-112.967 

********** 

-167.155 

167.155 

-263.589 

263.539 

********** 

0.000) 
0.000) 
**********) 
********** ) 
9.756) 

71 .633 
167 . 155 
86.330 
66 . 295 
9.756 

SL22T 

SL22C 

SLUT 

SL11C 

SL12S 


PLY FflILURE CRITERIA (NO temperature or 


MOISTURE STRESSES) 


PLY NO. 

= 3 

THETA = 90. 

■00 MATERIAL 

SYSTEM = SGLAHMHS AS— IMHS 



LOADS 

SLUT 

216.8321 

KSI 

SL11C 

166.5112 

KSI 

SL22T 

9.9151 

KSI 

SL22C 

23.1353 

KSI 

SL12S 

11.9513 

KSI 

FAIL. LOAD 
KSI 

MODE 

SCXXT MIN ( -5436.801 
SCXXC MIN ( 5436 .801 
SCYYT MIN C 86.330 
SCYYC MIN ( -86.330 
'SCXYS MIN (********** 

4175.063 
-4175.063 
-66 .295 
66 .295 
********** 

71.638 

-71.638 

112.967 

-112.967 

********** 

-167.155 

167.155 

-263.589 

263.589 

********** 

0.000) 
0.000) 
**********) 
**********) 
9.756 ) 

71.638 
167.155 
36 . 330 
66.295 
9.756 

SL22T 

SL22C 

SLUT 

SL11C 

SL12S 


. XTE; «****###***„ IMPLIES -« N0T APPLICABLE"- 


LAMINATE FAILURE STRESSES 


LAYUP — > 0 90 90 0 

FAILURE CRITERIA (NO TEMPERATURE OR MOISTURE STRESSES) 


PLY NO. 

= 4 

THETA = 0 . 

■00 MATERIAL 

SYSTEM = AS- -IMLS 

AS --IMLS 



LOADS 

SLUT 
22 2.7741 
KSI 

SL11C 
87 .6392 
KSI 

SL22T 

5.0065 

KSI 

SL22C 

15.0194 

KSI 

SL12S 

5.1261 

KSI 

FAIL. LOAD 
KSI 

MODE 

SCXXT MIN ( 
SCXXC MIN ( 
SCYYT MIN ( 
SCYYC MIN ( 
SCXYS MIN ( 

155.699 
-155.699 
-5076 . 309 
5076 . 309 
0.000 

-61.252 
61.252 
1997.017 
-1997 . 017 
0 .000 

251.062 
-251.062 
20 .504 
-20 . 504 
********** 

-753.187 

753.187 

-61.511 

61.511 

********** 

0.000 ) 
0.000 ) 
0.0 00 ) 
0.000 ) 
5.776) 

155.699 
61.252 
20.504 
61 .511 
5.776 

SL11T 

SLUC 

SL22T 

SL22C 

SL12S 
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ORIGINAL PAGE IS 
Item 20 OF POOR QUALITY 

S U PI M A R Y 

LAMINATE FAILURE STRESS ANALYSIS - (NO TEMPERATURE OR MOISTURE STRESSES) 

C BASED UPON FIRST PLY FAILURE) 


LOAD TYPE 

STRESS KSI 

FAILURE MODE 

PLY NO. 

THETA 

MATERIAL 

SYSTEM 

SCXXT 

71.638 

SL22T 

3 

90 . 0 

SGLAHMHS 

AS — IMHS 

scxxc 

61.252 

SL11C 

4 

0 . 0 

AS--IMLS 

AS — IMLS 

SCYYT 

20 . 50<+ 

SL22T 

1 

0 . 0 

AS--IMLS 

AS — IMI.S 

SCYYC 

61.511 

SL22C 

1 

0.0 

AS — IMLS 

AS--IMLS 

SCXYS 

5.776 

SL12S 

4 

0.0 

AS — IMLS 

AS — — IMLS 


LAMINATE FAILURE STRESS ANALYSIS - (NO TEMPERATURE OR MOISTURE STRESSES) 
(BASED UPON FIBER FAILURE) 


LOAD TYPE STRESS KSI FAILURE MODE PLY NO. THETA 

SCXXT 155.699 
SCXXC 61.252 
SCYYT 86.330 
SCYYC 66.295 
SCXYS ********* 


NOTE: IF THERE IS NO ANGLE PLY "SCXYS" BASED UPON FIBRE FAILURE IS NOT PREDICTED. 


SLUT A 0.0 AS--IMLS AS--IMLS 

SL11C 4 0.0 AS— IMLS AS-IMLS 

SL11T 2 90.0 S3LANI1US AS— IlltiS 

SL nc 2 9 0.0 SGLAHMHS A3 — IMHS 

N/A 
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Appendix C 

Resident Data Bank (FBMTDATA .B AN K) 


T300 

FP 3000 0.300E-03 0.640E-01 

FE 0.320E 08 0.200E 07 0.200E 00 0.250E 00 0.130E 07 

FT -0.550E-06 0.560E-05 0.580E 03 0.580E 02 0.170E 00 

FS 0.350E 06 0.300E 06 0.000 0.000 0.000 


0.700E 06 
0.000 


AS— 

FP 10000 0.300E-03 0.630E-01 

FE 0.310E 08 0.200E 07 0.200E 00 0.250E 00 0.200E 07 0.100E 07 

FT -0.550E-06 0.560E-05 0.580E 03 0.580E 02 0.170E 00 

FS 0.400E 06 0.400E 06 0.000 0.000 0.000 0.000 

SGLA 

FP 204 0.360E-03 0.900E-01 

FE 0.124E 08 0.124E 08 0.200E 00 0.200E 00 0.517E 07 0.517E 07 

FT 0.280E-05 0.280E-05 0.750E 01 0.750E 01 0.170E 00 

FS 0.360E 06 0.300E 06 0.360E 06 0.300E 06 0.180E 06 0.180E 06 

HttSF HIGH MODULUS SURFACE TREATED FIBER. 

FP 10000 0.300E-03 0.703E-01 

FE 0.550E 08 0.900E 06 0.200E 00 0.250E 00 0.110E 07 0.700E 06 

FT -0.550E-06 0.560E-05 0.580E 03 0.580E 02 0.170E 00 

FS 0.280E 06 0.200E 06 0.000 0.000 0.000 0.000 

OVER END OF FIBER PROPERTIES. 


IKLS INTERMEDIATE MODULUS LOW STRENGTH MATRIJ 


MP 

0.460E 

'01 





ME 

0.500E 

06 

0.410E 

00 

0.570E-04 


MT 

0.125E 

01 

0.250E 

00 



MS 

0.700E 

04 

0.210E 

05 

0.700E 04 

0.140E- 

MV 

0.225E 

00 

0.420E 

03 



IMHS INTERMEDIATE MODULUS HIGH STRENGTH 

MP 

0.440E- 

-01 





ME 

0.500E 

06 

0.350E 

00 

0.360E-04 


MT 

0.125E 

01 

0.250E 

00 



MS 

0.150E 

05 

0.350E 

05 

0.130E 05 

0.200E- 

MV 

0.225E 

00 

0.420E 

03 



HMHS HIGH 

MODULUS HIGH 

STRENGTH MATRIX. 

MP 

0.450E- 

■01 





ME 

0.750E 

06 

0.350E 

00 

0.400E-04 


MT 

0.125E 

01 

0.250E 

00 



MS 

0.200E 

05 

0.500E 

05 

0.150E 05 

0.200E-I 

MV 

0.225E 

00 

0.420E 

03 



OVER 

END 

OF 

MATRIX 

PROPERTIES . 
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